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ABSTRACT 


MICHAEL  L.  SHELLEY.  A  Model  of  Inhalation  Gas  and 
Vapor  Transport  in  the  Human  Lung.  (Under  the 
direction  of  DR.  ROBERT  L.  HARRIS) 


^The  transport  of  inhaled  gas  or  vapor  in  the  human  lung 
is  theoretically  modeled  to  determine  the  fraction  of 
inhaled  gas  or  vapor  absorbed  through  the  bronchial  walls  in 
the  conductive  zone  and  the  fraction  entering  alveolated 
space  and  undergoing  blood/gas  exchange.  An  iterative- 
analytical  approach  is  used,  rather  than  a  numerical  one,  to 
generate  a  large  amount  of  data  over  a  wide  range  of  gas  or 
vapor  solubilities  and  dif fusivities.  The  model  applies  to 
gases  or  vapors  of  low  solubility.  Simple  analytical 
expressions  are  fitted  to  model  results,  expressing 
fractional  bronchial  absorption  and  alveolar  uptake  in  terms 
of  gas  or  vapor  solubility  and  diffusivity.  Typical 
physiological  parameters  are  assumed  and  only  physical 
properties  of  the  gas  or  vapor  are  required  for  input. 
Methods  are  presented  which  allow  estimation  of  these  input 
values . 
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Chapter  1 


INTRODUCTION 

The  modeling  of  inhalation  exposure  to  gases  and  vapors 
has  received  much  attention  in  the  literature,  particularly 
in  the  anesthesiology  field.  Models  have  been  rigorously 
developed  for  alveolar  gas  exchange,  distribution  to  various 
tissue  types,  and  clearance  from  these  tissues  under  varying 
physiological  conditions.  Application  of  these 
mathematically  derived  relations  to  in-vivo  laboratory 
experiments  has  allowed  the  establishment  of  important 
physiological  constants  for  many  exposure  gases  (blood/gas 
and  tissue/blood  partition  coefficients,  tissue  clearance 
rates,  etc.).  These  data  provide  a  tool  for  predicting  body 
burdens,  half  life,  and  pharmacological  effects  of 
particular  gaseous  insults  at  various  inhalation  exposure 
concentrations.  In  large  part,  these  constants  are  based  on 
direct  measurement  of  blood  and  tissue  concentrations,  urine 
excretion,  ventilation,  etc.  as  well  as  certain  indirect 
measurements  such  as  analysis  of  end-exhalation  expired  air 
representing  alveolar  concentration.  Measurement  of 
alveolar  concentration  in  this  way  probably  provides  a 
reasonably  accurate  estimation  when  dealing  with  large 
exposure  concentrations  such  as  those  of  interest  in 


In  public  health,  one  is  more  often  concerned  with 
relatively  low  exposure  concentrations  of  many  types  of 


gases  with  a  wide  range  of  solubilities.  Under  these 
conditions,  distribution  of  gas  within  the  lung  must  be 
considered  in  a  different  light,  with  particular  attention 
to  the  significance  of  the  muco-ciliary  clearance  mechanism 
and  absorption  through  the  epithelial  wall  of  the  conductive 
zone  of  the  lung  prior  to  reaching  the  alveolar  region. 
Conceptually,  one  must  envision  a  mass  balance  which 
distributes  the  inhaled  mass  of  gas  into  the  known  pulmonary 
clearance  mechanisms  of  exhalation,  alveolar  gas  exchange, 
bronchial  wall  absorption,  and  mucus  discharge  to  the  GI 
tract.  Gas  exchanged  at  the  bronchial  epithelium  beneath 
the  mucus  blanket  should  ultimately  enter  the  same 
physiological  compartment  as  that  exchanged  in  the  alveoli 
but  with  a  different  rate  of  uptake  and  with  an  opportunity 
to  undergo  metabolism  in  lung  tissue  prior  to  entering  the 
blood  stream  for  further  distribution.  Sixteen  different 
enzymatic  systems  involved  in  metabolism  have  been  reported 
in  lung  tissue.  The  mixed  function  oxidase  systems  in  the 
lung  are  similar  to  those  in  the  liver.  Enzymatic  activity 
in  the  lung  is  small  relative  to  that  in  the  liver,  but 
considerable  metabolism  might  be  expected  due  to  higher 
exposure  to  inhalation  agents,  the  lung  acting  as  a  gateway 


An  understanding  of  the  distribution  partitioning 
described  above  is  essential  in  relating  exposure 
concentration  to  blood  concentration  as  well  as  estimating 
the  GI  insult  due  to  ciliary  clearance.  Yet,  such  a  model, 
using  this  mass  balance  approach,  to  date  has  not  been 
presented  in  the  literature.  The  complexities  of  gas 
exchange  back  and  forth  across  the  air/liquid  interface  on 
the  bronchial  wall,  together  with  the  uncertainties  of  the 
structure  and  composition  of  the  mucus  blanket  and  the  lack 
of  computer  methods  to  generate  data  consistent  with  lung 
morphology,  hindered  such  efforts  prior  to  the  early  50's 
when  existing  inhalation  models  were  developed. 

Construction  of  a  model  which  would  properly  distribute 
exposure  in  the  lung  relative  to  applicable  physical  and 
physiological  parameters  is  now  within  reach  and  would 
constitute  a  significant  contribution  to  the  study  of 
gaseous  inhalation  exposure. 

For  example,  consider  the  conceptual  model  in  figure 
1.1  which  is  of  the  type  often  used  by  systemic  distribution 
modelers.  In  this  model,  the  inhaled  mass  entering  the  dead 
space  (conductive  zone)  and  the  alveoli  (respiratory  zone) 
is  partitioned  into  the  two  compartments  of  blood  and 
exhalation.  The  mass  partitioned  into  the  blood  compartment 
is  based  on  complete  air/blood  equilibration  of  all  mass 
entering  the  alveolar  region  in  accordance  with  the  alveolar 
ventilation  rate  and  blood  perfusion  rate.  This  assumes 
that  the  mass  entering  the  alveolar  region  is  at  the  same 


F I GURE  1 . 1 


concentration  as  that  of  the  environment  (allowing  for  minor 

differences  due  to  body  temperature  and  humidity).  The 

equilibrated  blood  concentration  thus  arrived  at  is  the 

basis  of  subsequent  uptake  and  elimination  determination  in 

the  various  tissue  types.  (Note:  Ketty"^  presents  several 

approaches  to  this  type  of  model,  reviewing  previous  work  by 

30  10  26 

Zunts,  von  Schroder,  Widmark  ,  Haggard  ,  Teorell  , 

17  8 

Morales  and  Smith  ,  and  Copperman.  Fi serova -Be rgerova 

also  presents  a  fundamental  review  of  the  principles 

involved.)  For  the  case  of  low  concentration  exposure  to  an 

environmental  agent,  a  more  desirable  model  is  depicted  in 

figure  1.2  and  includes  a  more  detailed  partitioning  within 

the  lung.  Here,  rather  than  all  mass  passing  through  the 

conductive  zone  reaching  the  alveoli,  two  additional 

compartments  are  added:  bronchial  wall  absorption  to  lung 

tissue  and  mucus  discharge  at  the  top  to  the  GI  tract. 

Exposure  to  the  GI  tract  via  this  route  may  or  may  not  be 

significant.  Partitioning  to  lung  tissue  would  limit  the 

rate  of  transport  to  blood  and  provide  the  opportunity  for 

metabolism  prior  to  reaching  other  tissue  types. 

Partitioning  to  lung  tissue  or  discharge  to  GI  tract, 
however,  is  not  as  simple  as  assuming  complete  equilibration 
with  that  compartment.  One  must  describe  the  transport 
process  in  terms  of  the  gas  solubility  and  diffusivity 
through  the  medium.  If  such  a  description  can  be  achieved, 
the  behavior  of  the  gas  in  the  lung  could  be  predicted 
without  experimental  data  based  on  the  physical  constants  of 


the  gas.  The  relative  significance  of  exposure  via 
blood/gas  exchange,  bronchial  wall  absorption,  and  GI 
ingestion  by  mucus  discharge  could  then  be  analyzed  with 
confidence,  possessing  a  valuable  tool  in  estimating  the 
potential  toxicity  of  the  exposure. 

The  present  work  examines  the  fundamentals  involved  in 
creating  a  model  to  partition  inhaled  gas  or  vapor  into  the 
four  compartments  of  blood/gas  exchange,  bronchial 
epithelium  absorption,  mucus  discharge,  and  exhalation. 
First,  the  architecture  of  the  conductive  zone  of  the  human 
lung  is  examined  and  pathways  for  mass  transport  are 
considered.  Air  flow  patterns  in  light  of  this  structure 
are  studied.  Next,  the  theory  of  mass  transport  is 
presented  and  applied  to  human  lung  structure.  Differential 
equations  of  change  are  derived  for  various  transport 
regions  involved  in  lung  partitioning.  Next,  physical 
properties  of  inhaled  gases  and  vapors  are  discussed  in 
relation  to  their  impact  on  mass  transport  models.  Methods 
for  estimating  solubilities  and  dif fusivities  are  provided. 
Finally,  a  simplified  model  of  inhaled  mass  partitioning  in 
the  lung  is  presented,  using  an  iterative-analytical 
approach.  Results  of  the  model  over  a  range  of  exposure 
substances  allow  formulation  of  simple  expressions  to 
calculate  fractional  partitioning  into  compartments 
proceeding  from  the  lung. 


Chapter  2 

HUMAN  LUNG  STRUCTURE 


Any  mass  transport  model  must  make  precise  assumptions 
concerning  the  structure  of  the  transport  medium.  In  the 
human  lung,  this  can  be  quite  complicated  with  the  number  of 
bifurcating  bronchi  increasing  geome' rically  as  one  proceeds 
from  the  trachea  toward  the  alveolar  spaces.  The  length, 
diameter,  and  angle  of  bifurcation  of  the  bronchi  affect  air 
flow  patterns  and  transfer  of  mass  to  the  bronchial  walls, 
thus  affecting  absorption  rates.  These  dimensions  vary 
widely  from  bronchial  generation  to  generation  as  well  as 
from  individual  to  individual.  Once  absorbed  at  the 
bronchial  wall  surface,  further  mass  transport  of  a  gas  or 
vapor  into  various  fluid  and  tissue  layers  again  depends  on 
the  geometry  and  dimensions  of  these  regions  which  also  may 
have  wide  variations.  A  subset  of  structural  models  is 
required  which  will  characterize  the  dimensions  within  a 
particular  bronchial  generation  and  represent  an  average 


based  on  limited  experimental  data 


The  system  of  airway  branching  in  the  conductive  zone 
of  the  human  lung  serves  to  distribute  air  from  the  trachea 
to  the  300  million  alveoli.  The  airways  are  organized  in  a 
tree-like  system,  each  generation  branching  by  dichotomy 
(dividing  into  two  daughter  branches).  Figure  2.1  depicts 
this  bifurcating  system  and  illustrates  the  generation 
numbering  convention,  the  trachea  being  generation  'O'  and 
numbering  up  with  each  successive  bifurcation.  The 
conductive  zone  extends  to  approximately  generation  16,  the 
walls  of  which  are  relatively  smooth  consisting  of  varying 
degrees  of  cartilage  and  smooth  muscle.  Beyond  generation 
16  lie  the  transitory  and  respiratory  zones  where  bronchial 
walls  become  'spotted'  with  alveoli  becoming  more  and  more 
dense  until  the  walls  become  completely  alveolated.  Airways 
terminate  in  alveolar  sacs  at  approx im  itely  generation  23. 

Two  systems  of  dichotomy  are  depicted  in  figure  2.2: 
'regular'  and  'irregular'.  In  the  case  of  irregular 
dichotomy  the  angle  of  bifurcation,  diameter,  and  length  of 
the  two  daughter  bronchi  are  unequal,  whereas  with  regular 


dichotomy  the  bifurcation  is  completely  symmetric.  It  has 
20 

been  shown  that  the  human  lung  branches  by  irregular 

29 

dichotomy.  Weibel  has  presented  simple  expressions  for 


12 


for  a  regular  dichotomy  assumption.  Although  regular 
dichotomy  rarely  occurs  in  nature  and  is  less  representative 
of  reality,  it  is  more  simply  described  and  more  amenable  to 
computer  methods  employed  in  distribution  of  gases  in  the 
lung  than  in  irregular  dichotomy. 

Weibel's  regular  dichotomy  expressions  for  the  diameter 
and  length  of  airways  in  generations  0  thru  3  are 
exponential  functions  of  the  tracheal  diameter  and  length, 
respectively : 


°GEN  =  Vxp(-0.388  GEN)  (2.1) 

LGEN  =  L0exp(-°.92  GEN)  (2.2) 

where  DQ  =  tracheal  diameter  =  1.8  cm 

and  Lq  =  tracheal  length  =  12  cm 

Generations  4  thru  16  require  empirical  values  for  DQ  and 


°GEN  =  D' 0exp[-(° -293  -  0.0062  GEN ) GEN ]  (2.3) 

LGEN  =  L'oexP<-0*17  GEN)  (2.4) 


where  D' q  =  1.3  cm 

and  L  * q  =  2.5  cm 


Since  the  airways  branch  by  dichotomy,  the  total  number  of 
branches  within  a  generation  is  given  by 


N 


GEN 


2  N 


=  2GEN 


(GEN-1 ) 


(2.5) 
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Combinations  of  equations  2.1  thru  2.5  allow  complete 
description  of  the  gross  airway  dimensions  of  the  conductive 
zone.  Table  2.1  tabulates  some  key  dimensions  for  this 
reg ion . 

Air  Flow  Patterns 

The  structural  model  described  above  assumes  a  system 

of  ideally  bifurcating  tubes  of  ideal  shape  (right  circular 

cylinders).  Flow  patterns  within  such  systems  have  been 

18  23 

extensively  studied  and  discussed.  '  Flow  patterns,  of 
course,  depend  not  only  on  pathway  structure  but  on  local 
flow  velocity.  During  breathing  this  is  particularly 
complicating  due  to  the  oscillating  breathing  process.  As 
illustrated  in  figure  2.3  (curve  A),  the  flow  velocity  is 
near  zero  at  the  beginning  and  end  of  inhalation  or 
exhalation  and  maximum  near  the  middle.  Since,  during 
normal  breathing,  inhalation  is  generally  a  forced  process 
while  exhalation  is  passive,  the  curve,  in  reality,  may  be 
expected  to  be  slightly  skewed  in  varying  degrees  from 
inhalation  to  exhalation.  However,  for  purposes  of 
modeling,  one  may  consider  the  velocity  curve  to  be 
triangular  (curve  B)  with  the  maximum  velocity  at  mid 
inhalation  or  exhalation  dependent  on  the  tidal  volume  and 
respiratory  rate. 

Given  the  path  structure  and  local  velocity,  transport 
of  a  particular  gas  is  by  convection  (with  the  flow  of  air) , 


AIRWAY  DIMENSIONS  OF  CONDUCTIVE  ZONE 


(Regular  Dichotomy) 


NUMBER 


DI A( cm) 


LENGTH ( cm) 


TOTAL 
VOLUME  (cm) 


0 

1 

1.8 

12  .0 

30.50 

1 

2 

1.22 

4.76 

11.25 

2 

4 

0.83 

1  .90 

3.97 

3 

8 

0.56 

0.76 

1.52 

4 

16 

0.45 

1.27 

3.46 

5 

32 

0.35 

1.07 

3.30 

6 

64 

0.28 

0.90 

3.53 

7 

128 

0.23 

0.76 

3.85 

8 

256 

0 . 186 

0.64 

4.45 

9 

512 

0.154 

0.54 

5.17 

0 

1024 

0.130 

0.46 

6.2  1 

1 

2048 

0.109 

0.39 

7.56 

2 

4096 

0.09  5 

0.33 

9.82 

3 

819  2 

0.082 

0.27 

12.4  5 

4 

16384 

0.074 

0.23 

16.40 

5 

32768 

0.066 

0.20 

21.70 

6 

65536 

0.060 

0.165 

29.70 

TABLE 

2.1 
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molecular  diffusion,  and  absorption  at  the  surface  of  the 
tube  wall.  Molecular  diffusion  is  assumed  to  proceed  in 
accordance  with  Fick's  Law  (dependent  on  concentration 
gradients)  while  absorption  at  the  surface  depends  on 
concentration  gradients  and  uhe  solubility  of  the  gas  at  the 
prevailing  temperature.  Convection  (described  by  flow 
patterns)  influences  the  development  of  concentration 
gradients  and  thus  impacts  all  phases  of  transport. 

Within  the  convective  process,  flow  velocity  determines 
the  nature  of  the  velocity  profile  (i.e.  turbulent  or  non- 
turbulent).  Although  turbulent  conditions  can  occur  with 
some  extreme  physiological  conditions  in  the  first  three  or 
four  generations,  velocity  profiles  throughout  the  lung  are 
generally  considered  non-turbulent ,  taking  on  a  parabolic 
snape.  Figure  2.4  illustrates  some  aspects  of  this  type  of 
flow.  A  parabolic  velocity  profile  is  established  due  to 
frictional  resistance  at  the  walls  and  a  radial  transfer  of 
momentum.  There  is  near  zero  velocity  at  the  walls  and 
maximum  velocity  at  the  centerline  which  will  be  shown  in 
appendix  A  to  be  twice  that  of  the  average  velocity  in  any 
bronchial  generation.  Since  mass  is  convectively 
transported  at  different  rates  depending  on  radial  distance 
from  the  center,  concentration  gradients  controlling 
molecular  diffusion  are  thus  affected  to  varying  degrees. 
Additionally,  the  right  side  of  figure  2.4  illustrates  a 
further  complicating  factor.  Depending  on  the  angle  of 
bifurcation  and  assuming  that  the  velocity  profile  is 


TURBULENT  FLO 
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relatively  undisturbed  across  bifurcations,  an  axial 

asymmetry  may  be  expected  to  occur  in  the  daughter  branches 

as  the  velocity  profile  is  split  at  the  centerline  with  the 

maximum  velocity  in  the  daughter  branches  off-set  from 

center  toward  the  point  of  bifurcation.  Schroter  and 
23 

Sudlow  also  discuss  the  occurrences  of  vortices  in  the  air 

flow  pattern  at  bifurcations  which  appear  to  have  differing 

patterns  from  inhalation  to  exhalation. 

In  addition  to  the  above  aspects  of  convection,  studies 
22  27 

have  shown  '  that  longitudinal  dispersion  occurs, 
described  by  a  dispersion  coefficient  which  is  clearly 
greater  than  the  molecular  diffusion  coefficient  alone. 

This  effect  is  commonly  termed  'the  Taylor  effect'.  Taylor 
dispersion  is  simply  illustrated  in  figure  2.5.  Assume  an 
initial  plug  of  uniform  concentration  (hatched  area)  is 
introduced  at  the  entrance  of  a  cylindrical  tube.  In  the 
absence  of  molecular  diffusion,  a  parabolic  velocity  profile 
will  advance  the  plug  to  occupy  region  2  at  some  time  later 
with  zero  concentration  in  regions  1  and  3.  Without 
diffusion,  the  average  cross-sectional  concentration  will 
thus  reflect  significant  longitudinal  dispersion.  Molecular 

diffusion  will  further  spread  the  mass  from  region  2  towards 
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regions  1  and  3.  Taylor  presents  an  expression  describing 

the  average  concentration  at  a  distance  'X'  downstream  at 

time  't'  given  a  constant  uniform  concentration  CO 

introduced  at  X=0  with  a  centerline  velocity  v: 


t*V 


(2.6) 
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c/co  = 

0.5  +  0.5 

erf (0.5 

X1K"°*5t  °*5) 

(x1<=  0) 

c/co  = 

0.5  -  0.5 

erf (0.5 

X1K_0*5t"°*5) 

O 

/\ 

X 

where 

xx=  X 
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v2/( 19  2 

D) 

A  =  tube  radius 

D  =  molecular  diffusion  coefficient 


and  erf  =  the  error  function 


(2.7) 


This  formulation  assumes  that  the  radial  'smoothing  out'  of 
the  concentration  profile  due  to  molecular  diffusion  is 
rapid  compared  to  convective  effects.  It  also  assumes  no 
absorption  at  the  wall  which  limits  its  practical  use  in 
modeling  highly  soluble  mass  transport  in  the  human  lung. 
Taylor  also  cautions  that  this  expression  should  be  applied 
only  at  some  finite  time  and  distance  downstream  from  the 
tube  entrance  after  'stabilization'  of  the  assumed 
convective-diffusive  conditions.  Attempts  have  been  made  to 
circumvent  this  restriction  by  employing  time-dependent 
dispersion  coefficients.  For  a  detailed  discussion,  the 

9 

reader  is  referred  to  Gill  and  Sankarasubramanian  . 

However,  most  modeling  efforts  of  mass  transport  in  tubes 
assume  an  infinite  tube  length  ignoring  entrance  and  exit 
effects . 

An  excellent  summary  of  air  flow  patterns  in 
bifurcating  tubes  is  provided  by  Overton,  pages  93-102  of 
Miller  and  Menzel^^. 
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Structure  and  Composition  of  the  Bronchial  Wall 

Once  molecules  have  been  absorbed  by  the  bronchial 
wall,  further  absorption  is  influenced  not  only  by 
solubility  and  concentration  gradients  in  the  air  phase  but 
also  by  the  rate  of  transport  in  the  liquid  or  tissue  phase 
away  from  the  surface  and  deeper  into  the  wall.  Diffusion 
kinetics  are  active  here  as  well  and  the  structure  and 
dimensions  of  regions  beneath  the  surface  must  be  defined. 

The  surface  of  the  airways  is  bathed  with  a  fluid  known 
to  be  active  in  the  capture  and  clearance  of  particulate 
matter  in  the  conductive  zone.  Ciliated  cells  lining  the 
epithelial  wall  move  the  fluid  blanket  upward  toward  the 

tracheal  opening  where  it  is  expelled  to  the  GI  tract. 

2  8 

Wanner  has  published  a  detailed  review  of  the  structure 
and  function  of  this  muco-ciliary  clearance  mechanism.  The 
airway  epithelium  consists  of  several  types  of  ciliated  and 
non-ciliated  cells,  some  or  all  of  which  are  involved  in  the 
secretion  of  fluid  for  the  mucus  blanket.  The  relative 
number  of  ciliated  and  secretory  cells  generally  decrease 
from  the  trachea  downward  suggesting  a  systematic  variation 
from  generation  to  generation  of  fluid  thickness,  velocity, 
and  composition.  These  variations  have  not  been  clearly 
defined,  but  average  values  have  been  established  which  are 
perhaps  more  representative  of  the  upper  than  of  the  lower 


airways 


Figure  2.6  depicts  a  longitudinal  view  of  a  bronchial 
wall  segment.  The  cilia  (approximately  6  micro-meters  in 
length)  are  surrounded  by  a  watery  layer  of  fluid  just  above 
the  epithelial  wall.  Above  the  watery  layer  is  a  second 
fluid  layer  with  a  high  g lyco-protei n  content.  The  presence 
of  these  molecules  gives  this  layer  a  viscous  consistency. 
The  q lyco-protein  layer  rides  above  the  watery  layer  and  is 
propelled  toward  the  top  of  the  trachea  by  interaction  with 
the  tips  of  the  cilia.  The  direction  of  motion  of  the 
watery  layer  beneath  is  undefined  and  it  is  expected  to  be 
in  a  turbulent  state  due  to  the  rapid  and  sudden  motion  of 
the  cilia.  Particulate  matter  inhaled  into  the  conductive 
zone  which  is  captured  and  absorbed  by  the  g lyco-protein 
layer  would  be  expected  to  remain  captured  there  traveling 
with  the  layer  and  finally  expelled  to  the  GI  tract  perhaps 
hours  later.  Gaseous  molecules,  on  the  other  hand,  when 
captured  and  absorbed  into  the  g lyco-protein  layer  would 
continue  to  pass  into  and  out  of  various  fluid  or  tissue 
regions  (or  back  into  air)  depending  on  prevailing 
concentration  gradients  and  the  diffusivity  of  the  gas. 
Beneath  the  watery  layer,  the  epithelial  wall  itself 
consists  of  layers  of  surface,  intermediate,  and  basal  cells 
totalling  approximately  30  micro-meters  in  thickness  prior 
to  reaching  well-perfused  lung  tissue  where  the  gas 
concentration  would  be  in  equilibrium  with  blood  flow  to  the 


There  is  some  evidence  that  the  mucus  blanket  (glyco¬ 


protein  layer)  may  be  discontinuous  toward  the  distal  end  of 
the  respiratory  tract,  existing  as  patches  or  not  existing 
at  all.  Components  of  both  fluid  layers  may  be  secreted  and 
re-absorbed  into  the  epithelium  in  response  to  varying 
physiological  stimuli.  In  the  absence  of  concrete  data,  it 
is  best  for  modeling  purposes  to  assume  a  continuous  blanket 
of  uniform  thickness  within  each  bronchial  generation. 

Figure  2.7  summarizes  the  regions  involved  in  mass 
transport  of  inhaled  gases  within  a  single  bronchial 
segment.  A  concentration  profile  is  illustrated  which,  in 
this  case,  depicts  a  gradient  from  the  air  phase  toward  lung 
tissue.  Molecules  are  convectively  transported  in  and  out 
of  the  lung  with  varying  velocity  in  the  air  (lumen)  while 
experiencing  radial  diffusion  as  well.  Once  absorbed  by  the 
g ly co-prote i n  layer  predominant  transport  is  by  molecular 
diffusion  with  a  distinct  concentration  profile  likely 
developing  within  this  viscous  layer.  A  slight  convective 
component  would  also  be  involved  here  as  the  layer  is  cilia- 
propelled  intact  at  a  rate  of  approximately  3  to  10  mm/min. 
Rapid  turbulent  mixing  in  the  watery  layer  would  predominate 
over  any  diffusive  effects,  resulting  in  a  uniform 
concentration  profile.  Vvithin  the  epithelium,  again  a 
distinct  profile  develops  depending  on  the  relatively 
uniform  tissue  concentration  beyond. 

The  suggested  structural  model  above  should  allow  a 
practical  description  of  transport  in  the  conductive  zone  of 
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the  human  lung  given  knowledge  of  the  contaminant  gas 
diffusivity  in  the  various  regions  and  its  solubility  which 
would  define  the  partitioning  of  gas  between  air  and  liquid 
phases  at  the  air-liquid  interface.  Development  of 
transport  expressions  and  models  in  subsequent  chapters  will 
refer  to  the  structural  model  depicted  in  figure  2.7. 


Chapter  3 
MASS  TRANSPORT 


As  demonstrated  in  chapter  2,  the  airways  of  the  human 
lung  are  arranged  in  a  complicated  tree-like  expansion, 
bifurcating  into  thousands  of  daughter  branches  of  varying 
dimensions.  Simplifying  assumptions  allow  consideration  of 
a  'regular  dichotomy'  model  in  which  all  branches  of  a 
particular  generation  have  identical  dimensions.  In  this 
way  a  relatively  simple  algorithm  may  be  employed  to  model 
the  transport  of  a  gas  or  vapor  simply  relating  the 
concentration  leaving  a  generation  to  the  concentration 
entering  an  adjacent  generation.  Prior  to  this  step, 
however,  transport  must  be  described  in  detail  on  the  level 
of  a  single  bronchial  segment,  in  accordance  with  a 
reasonable  structural  model  like  that  depicted  in  figure 
2.7.  This  chapter  will  present  equations  of  change  which 
describe  the  transport  phenomenon  within  each  transport 
reg ion . 

To  describe  transport  in  the  air  phase  (lumen  of  the 
bronchiole),  consider  the  elemental  volume  shown  in  figure 
3.1.  We  will  take  the  elemental  volume  to  be  a  ring  of 
height  Az,  depth  Ar,  and  width  w  =  2  77"r.  If  we  assume 
that  there  is  no  net  mass  flux  in  the  w  direction,  then  mass 


enters  and  leaves  the  elemental  volume  in  the  r  direction  by 
molecular  diffusion  (given  non-turbulent  air  flow)  and  in 
the  z  direction  by  molecular  diffusion  and  convection. 

Using  the  notation  'Nz|z'  to  denote  mass  flux  in  the  z 
direction  evaluated  at  the  point  'z',  we  can  write  the 
following  mass  balance  for  the  elemental  volume: 


[Nz  2  TT  r  Ar]  |z  -  [Nz  2  77  r  Ar]  |  z+  az  +  [Nr  2  77  r  az]  |r 
[Nr  2  77rAz]jr+Ar  =  d  C/  d  t  [27Tr  Ar  Az]  (3.1) 

The  right  hand  side  of  the  equation  merely  conserves  mass. 

If  the  sum  of  the  mass  entering  and  leaving  the  element  does 
not  equal  zero,  the  difference  must  reflect  a  change  in 
total  mass  over  the  period;  thus,  a  change  in  concentration 
(C)  over  time  times  the  volume  of  the  element.  Dividing 
equation  3.1  by  2  7Tr  Ar  az  ,  we  have 


Nz  I  z  -  Nz I z+  AZ 


+  [Nr  rl  |  r  -  [Nr  r]  I  r+  at  = 


r  Ar 


(3.2) 


The  terms  on  the  left  we  recognize  as  the  definition  of  the 
differential  if  we  take  the  limit  as  az  and  Ar  approach 
zero.  Thus  we  have 


-  ^  N  z  -  1  ^  ( Nr  r)  =  &C 


If  we  neglect  molecular  diffusion  in  the  z  direction  as 
insignificant  compared  to  convection,  we  may  describe  the 
mass  flux  as 


Nz 

=  v  C 

where 

V 

=  air 

flow  velocity 

(3, 

.4) 

and,  by 

Fick ' s 

Law 

Nr 

=  -DA 

bC/  b  r 

(3, 

.5) 

where 

DA 

=  the 

gas  diffusion  coefficient 

in  air. 

Substituting  equations  3.4  and  3.5  into  equation  3.3 

-  v  d  C  +  DA  b  [r  be/  b  r]  =  b  C 

dz  r  b  r  a  t  (3.6) 

or 

-v  be/  d  z  +  DA  d2C/  dr2  +  DA  bC/r  b  r  =  b  C/  b  t 

(3.7) 


Now  we  have  a  differential  transport  equation  in  terms  of 


known  parameters,  except  that  in  non  turbulent  flow,  the 


is  parabolic,  varying  with  r.  For  the  case 
through  a  circular  tube  this  profile  can  be 
of  a  radial  momentum  balance  and  can  be 


where  R  =  the  tube  radius  ( RG  in  figure  3.1) 

and  vAVG  =  *:'rie  v°lumetric  a^r  flow  rate  divided  by  the 

'  tube  cross-sectional  area 

(For  a  complete  derivation,  see  Appendix  A).  Finally, 
substituting  in  equation  3.7,  we  have 

-2  vAVG[l  -  (r/RG)2]  dc/  dz  +  DA  d2C/  dr2 

+  DA  dc/r  dr  =  dc/  bt  (3.9) 

which  describes  mass  transport  in  the  air  phase.  For  the 
case  of  turbulent  air  flow  (rare  in  the  human  lung),  the 
equation  becomes  much  more  complicated  and  must  be  divided 
into  two  equations,  one  applicable  to  the  turbulent  core  and 
one  to  the  region  near  the  wall.  For  a  derivation  of  the 
turbulent  case,  the  interested  reader  is  referred  to 
Appendix  B. 

For  transport  in  the  viscous  glyco-protein  region  [RG  < 
r  <  RW]  we  again  refer  to  an  elemental  volume  (figure  3.2). 
We  begin  with  the  mass  balance  of  equation  3.1  reducing  to 
the  differential  form  of  equation  3.3  as  in  the  case  of  the 
air  phase.  Ciliary  action  in  this  region  propels  the  mucus 
blanket  upward  (in  the  z  direction)  at  a  rate  less  than  0.5 
cm  per  minute.  This  is  a  slow  clearance  process  and  may  be 
considered  distinct  from  the  rapid  transport  process  of  air 
convection  and  molecular  diffusion.  Thus,  the  concentration 
in  the  glyco-protein  region  is  predominantly  influenced  by 
the  radial  concentration  gradient  due  to  diffusion,  and 
ciliary  convection  in  the  z  direction  is  neglected.  Thus, 


and 


Nz  =  -DG  d  C/  d  z 
Nr  =  -DG  dc/  dr 


(3.10) 


In  the  water  layer  [RW  <  r  <  RE] ,  rapid  turbulent 
mixing  distributes  any  entering  mass  uniformly  across  the 
region  instantaneously  (given  the  small  distance  across). 
Since  the  height  of  the  region  in  the  z  direction  is  large 
compared  to  (RE  -  RW) ,  diffusion  would  remain  a  factor  in 
that  dimension.  Thus,  we  consider  an  elemental  ring  of 
height  z  and  width  (RE  -  RW)  [figure  3.2]  and  write  the 
mass  balance  as  follows: 


[Nz  Az]  |z  -  [Nz  Az]  |  z+  az  +  [Nr  Ar]  |  RW 

-  [Nr  Ar)  |RE  =  dc/ dt  V  (3.13) 

where  Nz  =  -DW  dC/  dz  (3.14) 

Nr  |  RW  =  -DG  <)C/  <)r|RW  (3.15) 

Nr  |  RE  =  -DE  dC/drjRE  (3.16) 

Az  =  7 T  (RE2  -  RW2)  (3.17) 

Ar|RW  =  2  TT  RW  iz  (3.18) 

Ar  |  RE  =  2  77”  RE  az  (3.19) 

V  =  TT  (RE2  -  RW2)  Az  (3.20) 


DE  =  diffusion  coefficient  in  the  epithelium 


and 


DW  =  diffusion  coefficient  in  water 


Note  that  the  flux  in  and  out  in  the  radial  direction  is 
coupled  to  transport  in  the  adjacent  regions  and  assumes 
conservation  of  mass  across  the  boundaries.  Forming  the 
differential  in  the  z  direction  and  substituting  in  equation 
3.13,  we  have 

DW  d2C/  dz2  -  DG  [2  RW/(RE2  -  RW2  )  ]  b  C/  b  r  |  RW 
DE  [2  RE/(RE2  -  RW2  )  ]  b  C/<$r|RE  =  bC/bt 

(3.21) 

describing  transport  in  the  water  layer. 

The  equation  of  change  for  transport  in  the  epithelium 
is  directly  analogous  to  that  for  the  g lyco-protein  region 
and  may  be  written  immediately  by  inspection  of  equation 
3.12 

DE  b2C/  dz2  +  DE  b2 C/  b  r2 

+  DE  b  C/r  b  r  =  bC/  bt  (3.22) 

To  summarize,  the  final  equations  of  transport  for  each 
region  are  re-written,  representing  a  system  of  simultaneous 
partial  differential  equations  which  must  be  solved: 

Air: 

-2  VAVG11  "  (r/RG)2  ^c/  +  DA  <*2C/ 

+  DA  bc/r  be  =  be/  bt 


(3.9) 


Glyco-protein: 

DG  J2C/X  b  z2  +  DG  b2C/X  b  r2 

+  DG  bc/Xrbr  =  dC/X  e>  t  (3.23) 

Water : 

DW  b2C/X  b  Z2  -  DG  [2  RW/(RE2  -  RW2 )  ]  <3c/A  <3  r  |  RW 

DE  [2  RE/ (RE2  -  RW2  )  ]  b<Z/X  ^  r |  RE  =  C/A  <3  t 

(3.24) 

Epithelium: 

DE  b2C/X  b  z2  +  DE  <)2C/A  Cl  r2 

+  DE  <3  C/r  A  <3  r  =  b  C/a  <3  t  (3.25) 

Here,  A  is  an  air/liquid  partition  coefficient  which  allows 
C  to  be  expressed  in  terms  of  air  concentration  throughout. 
The  partition  coefficient,  like  the  diffusion  coefficient, 
is  a  characteristic  of  the  transported  gas  and  is  related  to 
its  solubility.  The  above  equations  assume  that  the  gas  is 
partitioned  across  the  air/liquid  interface  by  the  partition 
coefficient  (i.e.  air/liquid  equilibrium  exists  at  the 
surface).  Solution  of  the  above  equations  requires 
specification  of  boundary  conditions  which  include 
description  of  the  concentration  entering  the  bronchiole, 
zero  radial  flux  at  the  centerline,  conservation  of  mass 
across  region  boundaries,  and  constant  tissue  concentration. 
The  complexity  of  the  system  (even  under  simplifying 
assumptions  under  which  it  was  derived)  precludes  the 
development  of  a  single  analytical  expression  representing  a 
solution.  Numerical  analysis  would  generally  be  employed  to 
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study  the  transport  of  one  or  two  inhalation  gases.  A 
numerical  approach  is  not  trivial  and  would  consume  a  great 
deal  of  computer  resources.  Alternatively,  to  gather  a 
large  amount  of  transport  data  over  a  broad  range  of 
inhalation  gases  and  vapors,  further  simplifying  assumptions 
can  be  made  allowing  the  use  of  analytical  expressions  to 
describe  portions  of  the  transport  problem  and  approximating 
overall  results  in  an  iterative-analytical  approach.  Such  a 
model  is  described  in  chapter  5. 
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PHYSICAL  PROPERTIES  OF  INHALED 
GASES  AND  VAPORS 


Regardless  of  the  method  used  to  describe  mass 
transport  in  the  conductive  zone  of  the  human  lung,  any 
model  must  depend  on  accurate  values  of  inhaled  gas  or  vapor 
physical  characteristics:  namely,  its  solubility  and 
diffusivity.  Ideally,  one  would  hope  to  be  able  to 
reference  experimentally  derived  values  for  a  particular 
species  in  air  and  water.  We  assume  here  that  solubility 
and  diffusivity  in  the  mucus  layers  and  tissue  is  reasonably 
equivalent  to  that  in  water,  water  being  the  predominant 
constituent  of  composition.  Attempts  to  estimate  more 
accurate  values  based  on  an  analysis  of  all  possible 
chemical  interactions  between  solute  and  small  fractional 
components  of  the  solvent  would  be  difficult  to  justify  in 
light  of  the  benefit.  Because  such  values  are  not  readily 
available  in  the  literature,  reasonable  estimates  must  be 
used  unless  one  has  the  resources  available  to  gather  this 
data  experimentally.  This  chapter  presents  suggested 
methods  for  estimating  contaminant  gas  solubility  and 
diffusivity  with  specific  examples. 


In  mass  transport  models  describing  transport  across 
air/liquid  interfaces,  we  assume  that  the  concentration  at 
the  interface  just  to  the  liquid  side  is  in  equilibrium  with 
the  interface  concentration  just  to  the  air  side  and  that 
the  ratio  of  these  two  concentrations  is  constant  at  a  given 
temperature  and  pressure  (within  practical  inhalation 
concentration  ranges).  Thus  we  consider  a  partition 
coefficient  (oc)  which  converts,  for  example,  milligrams  per 
liter  of  air  (CA)  at  the  interface  to  milligrams  per  liter 
of  water  (CW)  at  the  interface: 

CW  =  oc  CA  (4.1) 

Definition  of  this  partition  coefficient  is  less  complicated 
in  the  lung  since  lung  temperature  and  pressure  do  not  vary 
enough  to  significantly  affect  it.  Thus,  oc  (for  an 
air/water  interface)  may  be  considered  to  depend  solely  on 
the  particular  gas  or  vapor  in  question. 

The  partition  coefficient  is,  of  course,  an  indicator 
of  the  gas's  solubility.  Estimation  of  solubilities  is  not 
well  developed  due  to  the  difficulty  in  describing  all 
possible  complicating  interactions  of  solute  and  solvent 
molecules.  In  the  ideal  case,  one  would  expect  the 
air/liquid  equilibrium  to  follow  Raoult's  Law: 

Yi  P  =  Xi  Pvp  (4.2) 

where  Yi  =  mole  fraction  of  solute  in  air 


Xi  =  mole  fraction  of  solute  in  liquid 

P  =  total  pressure  of  system 

and  Pvp  =  saturation  vapor  pressure  of  solute 

at  the  system  temperature 

This  law,  illustrated  by  the  solid  line  in  figure  4.1, 
simply  states  that  the  partial  pressure  exerted  by  solute  in 
air  approaches  linearly  the  saturation  vapor  pressure  of  the 
solute  as  the  mole  fraction  in  liquid  approaches  1.0  (pure 
solute).  Molecular  interactions  between  solute  and  solvent 
cause  the  equilibrium  function  to  depart  from  the  ideal. 

This  is  particularly  true  in  the  case  of  a  polar  solvent 
like  water  in  which  case  the  function  may  appear  as  that 
represented  by  the  dashed  line  in  figure  4.1.  Thus,  we 
require  a  correction  factor  /,  know  as  the  activity 
coefficient  which  will  adjust  the  slope  of  the  curve  from 
the  ideal  case  for  low  concentrations  encountered  in  the 
lung,  such  that 

Yi  P  =  /  Xi  Pvp  (4.3) 

21 

Reid,  Prausnitz,  and  Sherwood  review  several  methods  of 
estimating  activity  coefficients  for  vapor-liquid 
equilibria.  Combinations  of  empirical  equations  and 
equation  coefficients  are  offered  for  classes  of  solutes  and 
solvents  in  table  8-17  of  their  work.  Table  8-21  also 
provides  parameters  for  a  more  complicated  series  of 
calculations  in  a  functional  group  additive  contribution 
method.  Both  methods  are  semi-empirical:  the  first  assumes 


that  compounds  of  the  same  class  will  behave  similarly  while 
the  second  assumes  that  the  contribution  of  functional 
groups  is  additive  regardless  of  their  combination  with 
other  groups.  These  methods  apply  only  to  vapor-liquid 
solubility  and  not  to  solubility  of  substances  which  are 
gaseous  at  standard  conditions. 

13 

For  the  case  of  gases,  Hayduk  provides  direct 
empirical  values  for  some  g^ses.  His  values,  termed 
hydrogen-bonding  factors  ( cr  )  ,  effectively  function  as  the 
inverse  of  the  activity  coefficient.  That  is, 

Yi  P  =  Xi  Pvp/<7“  (4.4) 

where  a~  (H-bonding  factor)  =  the  ratio  of  Xi  to  the 

predicted  mole  fraction  in 
water  according  to  Raoults' 
Law. 

Appropriate  tables,  adapted  from  Reid,  Prausnitz,  and 
21  13 

Sherwood  and  from  Hayduk  are  reproduced  m  Appendix  G. 

E x amp le :  Determine  the  partition  coefficient  (  oc  )  for 
1 ,1 ,1-trichloroethane  and  propane  at  an  air/water  interface 
at  37  degrees  C  and  1  atm  pressure. 

|  If  CW  =  mg/1  in  water  and  CA  =  mg/1  in  air, 

|  then  we  seek  oc  satisfying: 

|  CW  =  oc  CA 

i 

Yi  =  CA/(1000  mol.wt.)  [gm-moles  per  liter] 

|  *  25.4  [liters  per  gm-mole  at  37  deg  C] 

|  =0.0254  CA/mol.wt. 
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Since  1  liter  of  water  is  55.49  gm-moles, 

CW  =  Xi  (55.49)  mol.wt.  (1000) 

=  55490  (mol.wt.)  Xi  (mg/1) 

and 

Xi  =  CW/(55490  mol.wt.) 

For  vapors, 

/  =  Yi/( Xi  Pvp )  =  0.0254 )CA/mol.wt. ) 

Pvp  CW/(55490  mol.wt.) 

=  1409  CA/(CW  Pvp) 

CW  =  [  124  3/(  if  Pvp)]  CA 
Thus,  oC  =  1409/(  Pvp) 

For  gases, 

<7-  =  Xi  Pvp/Yi  =  CW  Pvp/(14  09  CA)  [inverse  of/1] 

CW  =  [  14  09  O-  /Pvp]  CA 

Thus,  =  14  09  r  /Pvp 

For  tr ichloroethane  (a  mono-alkyl  chloride),  table  8- 

2  1 

17  of  Reid,  Prausnitz,  and  Sherwood  allows 

calculation  of  the  activity  coefficient  at  20  degrees 

C  as  380.0.  This  is  acceptable  as  a  good  estimate  for 

37  degrees  C  as  well  since  the  difference  would  be 

1  3 

insignif icant  as  shown  by  Hayduk  .  The  saturation 
vapor  pressure  (0.33  atm)  is  readily  available 


I  handbook  data  (interpolated  for  37  degrees).  Thus, 
|for  1 ,1 ,1-trichloroethane 

I 

c<  =  1409/(380.0  x  0.33)  =  11.2 

I 

| For  propane,  the  hydrogen  bonding  factor  is 

|read  directly  from  Hayduk  (table  2)  as  0.0326  and, 
|from  handbook  data,  Pvp  =  14.8  atm  (interpolated  to 
|  37  deg  .  )  .  Thus 

|  <5C  =  1409  (0.0326/14.8)  =  3.1 

Air  Diffusion  Coef f ic ient 

Methods  for  estimation  of  diffusion  coefficients  are 

slightly  more  developed  than  for  solubility.  For  air, 

2 1 

several  methods  can  be  found  in  the  literature  ,  all  of 

.  19 

which  recognize  the  same  fundamental  components.  Panwitz 

presents  the  following  equation: 

DFA  =  0.941077  T1*5/[P  W  (0.5  dft+  0.5  dg)2] 

*  (1/Ma  +  1/Mb)0*5  (4.5) 

2 

where  DFA  =  diffusion  coeff.  in  air  (cm  /sec) 

T  =  310  degrees  K  (body  temperature) 

Ma  =  28.96  (mol.wt.  of  air) 

Mn  =  mol.wt.  of  gas 

D 

d^  =  3.617  (dia  of  avg  air  molecule,  angstr) 

P  =  1013.3  (pressure  in  mbar) 

1  /3 

dn  =  dia  of  gas  molecule  =  1.18  VA 

D 
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VA  =  molar  volume  at  boiling  point  =  0.285 
VC  =  critical  volume  of  gas 
and  W  =  collision  integral. 

2  1 

Reid,  Prausnitz,  and  Sherwood  presents  the  following 
expression,  by  Neufeld,  for  the  collision  integral: 

W  =  A/SB  +  C/exp(  DS)  +  E/exp(  FS)  +  G/exp(HS)  (4.6) 

where  S  =  KT/f 

K  =  Boltzman's  constant 
£  =  characteristic  molecular  energy 

and 

A  =  1.06036;  B  =  0.15610; 

C  =  0.19300;  D  =  0.47635; 

E  =  1.03587;  F  =  1.52996; 

G  =  1.76464;  H  =  3.89411. 

S  can  also  be  estimated  in  terms  of  the  gas  boiling  point 
(BP)  as 

S  =  T/ (111.55  BP) 1/2  (4.7) 

Thus,  diffusivity  in  air  can  be  calculated  from  the 

molecular  weight  and  boiling  point  (handbook  data)  and  the 

critical  volume,  values  of  which  are  tabulated  for  a  wide 

2  1 

range  of  compounds  in  Reid,  Prausnitz,  and  Sherwood 

W a  ter  Dif  fusion  Coef  f ic ient 

The  best  and  most  widely  used  method  for  estimating 

31 

diffusivity  in  water  ( DFW )  is  the  Wilke-Chang  method 


.,TV 
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given  by 

DFW  =  7.4  X  lO"8  ($  M)1//2T/(  P?  VA°*6)  (4  'M 

where  M  =  18.02  (mol.wt.  of  water) 

T  =  3 10  degrees  K  (body  temperature) 

y  =  0.6915  cp  (viscosity  of  water) 

VA  =  gas  molar  volume  (estimated  from  the  critical 
volume  as  in  the  previous  section) 

and  $>  =  2.6  (association  factor  for  water) 

Example :  Determine  the  diffusion  coefficient  in  air  and 

water  for  benzene  (MW  =  78)  at  37  deg  C. 

| From  handbook  data,  the  boiling  point  (BP)  for  benzene 

2 1 

| is  80  deg  C  and  from  Reid,  Prausnitz,  and  Sherwood 
(the  critical  volume  (VC)  is  259  cc/gm-mole. 

I 

For  air: 

I 

| By  eq 1 n  4.7, 

|  S  =  310/[ 111.55(353) ] 1/2 
|  =1.56 

I 

and ,  by  eq 1 n  4.6, 

I 

W  =  1.06036/1.560.15610  + 

|  0.193/exp[0. 47635(1. 56)]  + 

1  .0  3587/exp[l. 52996(1. 56)]  + 

|  1.76464/exp[3. 89411(1. 56)] 

|  =  1.18 

VA  =  0.285(259)1,048  =  96.38 

I 

I  dB  =  1.18 (96.38 ) 1/3  =  5.41 


flat*-- 


DFA  =  0.941077(310) 1 *5/ 

[1013. 3(1. 18)  [0.5(3.617)  +  0.5(5.41)]  * 

'  0 
( 1/28.96  +  1/78] 

!  2 

=  4.6E-2  cm  /sec 

I 

|  For  water: 

|  By  eq ’ n  4.8, 

|  DFW  =  7.4E-8  [2.6(18.02)] 1/2310/ 

|  [0. 6915(96. 38  ) 0  *  6] 

2 

|  =  1.46E-5  cm  /sec 

In  summary,  mass  transport  modeling  in  the  conductive 
zone  of  the  lung  requires  values  of  air/liquid  partition 
coefficients  (  oC  or  PN),  air  diffusivity  (DFA),  and  water 
diffusivity  (DFW).  In  the  absence  of  reliable  experimental 
data,  these  values  can  be  estimated  by  existing  methods 
using  readily  available  handbook  data  as  illustrated  in 
figure  4.2.  Table  4.1  lists  solubility  and  diffusivity  data 
calculated  by  methods  described  in  this  chapter  for  some  of 
the  more  common  gases  and  vapors. 


Table  4.1 


WATER  SOLUBILITY  AND  DIFFUSIVITY  DATA  FOR  SOME 
COMMON  GASES  AND  VAPORS 


PARTITION 

cm^/sec 

GAS/VAPOR 

COEFF.  (PN  ) 

DIFF. COEFF 

Acetone 

348.0 

1.67E-5 

Acetylene 

0.81 

2.5E-5 

Ammonia 

237.0 

3.26E-5 

Benzene 

4.6E-6 

1 . 46E-5 

Butane 

17.2 

1.48E-5 

Carbon  Tetrachloride 

51.1 

1 .41 E-5 

Chloroform 

32.2 

1.54E-5 

Dichlore thane 

3.04 

- 

Dichloroethy lene 

5.32 

- 

Ethane 

0.034 

2.08E-5 

Heptane 

0.08 

1.06E-5 

Hexane 

0.11 

1.17E-5 

1 , 1 , 1 -Trie hloroe thane 

11.2 

1.35E-5 

Methylene  Chloride 

16.6 

1 .76E-5 

Methyl  Ethyl  Ketone 

378.0 

1.44E-5 

Nitrous  Oxide 

0.48 

2.71 E-5 

Octane 

0.05 

9.77E-6 

Pentane 

0.18 

1 .33E-5 

Pe  rc  hloroe  thy  le  ne 

74.0 

1.36E-5 

Propane 

3.1 

1 .71E-5 

Sulfur  Dioxide 

24 .0 

2.35E-5 

Toluene 

4 . 0E-6 

1 . 29E-5 

Trichloroethylene 

18.5 

1.47E-5 

Vinyl  Chloride 

1.6 

1.9E-5 

Xylene 

3.2E-6 

1.17E-5 

Chapter  5 

A  LUNG  MASS  TRANSPORT  MODEL 


As  discussed  in  chapter  3,  an  iterative-analytical 
model  is  desired  which  would  track  the  transport  of  an 
inhaled  gas  or  vapor  in  the  conductive  zone  of  the  lung  over 
a  sufficiently  long  exposure  time  to  determine  routes  of 
uptake.  Specifically,  over  the  course  of  an  exposure,  we 
wish  to  partition  a  volume  of  gas  or  vapor  into  fraction 
entering  system  through  alveolar  blood/gas  exchange, 
fraction  entering  system  through  bronchial  absorption  in  the 
conductive  zone,  fraction  expelled  to  the  GI  tract,  and 
fraction  exhaled.  If  a  large  amount  of  such  data  is 
generated  over  a  range  of  gases  and  vapors  of  varying 
physical  characteristics,  then  these  compartmental  uptakes 
can  be  characterized  as  a  function  of  solubility  and 
diffusivity,  assuming  that  physiological  conditions  remain 
constant . 

The  model  must  start  with  an  exposure  concentration 
entering  the  trachea,  distribute  that  mass  within  the 
trachea,  and  determine  the  concentration  leaving  and 
entering  the  next  generation.  A  concentration  grid  must  be 
established  within  each  region  of  transport  (figure  2.7) 
which  will  control  subsequent  mass  flux  as  inhalation 


continues.  The  time  of  inhalation  must  be  tracked  to 
determine  flow  velocity  and  rate  of  convection.  During 
exhalation  the  input  concentration  to  the  conductive  zone  is 
the  alveolar  concentration  which  is  based  on  air/blood 
equilibration  of  the  mass  entering  the  alveoli  during 
inhalation.  One  must  also  consider  that  the  tissue 
concentration  will  eventually  be  influenced  by  rising  blood 
concentrations  from  alveolar  exchange. 

We  begin  by  considering  the  transport  phenomenon  on  the 
level  of  a  single  bronchiole  as  shown  in  figure  2.7. 
Transport  in  the  air  phase  is  the  most  difficult  to  describe 
analytically.  Air  enters  at  one  end  with  an  existing 
concentration  profile.  As  the  model  is  moved  toward  over  a 
time  increment,  the  profile  is  influenced  by  laminar  flow 
convection  and  radial  diffusion  with  flux  through  the  wall 
while  experiencing  a  continually  changing  boundary  condition 
at  the  air/liquid  interface  as  it  moves  through  the 
bronchiole.  This  is  a  complex  analytical  problem,  and  we 
must  resort  to  considering  a  very  small  time  increment  over 
which  the  air  advances  only  a  small  distance  such  that  the 
condition  of  the  liquid  surface  is  considered  constant. 

With  this  approach,  we  can  consider  an  analytical  expression 
describing  the  profile  change  for  radial  diffusion  only, 
with  constant  boundary  condition,  and  then  adjust  that 
profile  according  to  the  convective  velocity  distribution  as 
shown  in  figure  5.1.  The  shape  of  the  profile  moved  forward 
would  depend  on  the  profiles  in  upstream  incremental  volumes 


and  the  fraction  of  those  volumes  moved  foward  at  each 

radial  point  by  the  parabolic  velocity  distribution.  The 

radial  diffusion  expression  would  be  analogous  to  other 

transport  problems  in  the  literature  such  as  heat  flow  in  a 

14 

cylindrical  rod.  Jaegar  presents  a  heat  conduction 
expression  defining  a  profile  as  a  function  of  distance  from 
centerline  [r],  and  initial  profile  [ f ( r ) ] .  In  terms  of 
mass  transport, 

C  =  Cl  [  1  -  2  J  0  (  r  <*»  )  e  x  p  (  -  DA  oc„2t)]  + 

RG  n=l  oc„  JO  (RG  °c„  ) 

2/RG2  JO  (r  )/Jl2(RG  )  exp(-DA  <*»2t)  * 

n=l 

r'f(r')  JO  (r 1  ocn  )dr'  (5.1) 

where  Cl  =  constant  surface  concentration 
t  =  time  increment 

J0,J1  =  zero  and  first  order  Bessel  functions  of  the 
first  kind 

and  the  oC„'s  are  the  real,  positive  roots  of  JO. 

Use  of  the  Bessel  functions  is  necessitated  by  the 
cylindrical  geometry  of  the  problem.  Unfortunately, 

Jaeger's  expression  (formulated  for  heat  conduction 
problems)  does  not  provide  satisfactory  results  for  many 
values  of  DA,  t,  and  RG  which  are  encountered  in  the  lung 
mass  transport  problem. 

However,  if  we  consider  the  case  of  a  gas  or  vapor 
which  is  not  highly  soluble  such  that  wall  absorption  is 


small,  then,  over  the  length  of  a  single  bronchiole,  wall 
flux  would  not  significantly  influence  the  concentration 
profile  in  air.  If  we  further  assume  that  mixing  at 
bifurcations  is  sufficient  to  uniformly  distribute  the  mass 
over  the  cross-section  at  entry,  then  a  relatively  uniform 
air  concentration  profile  may  be  assumed  within  each 
bronchiole  as  shown  in  figure  5.2.  In  this  way,  an  entire 
bronchiole  may  be  treated  as  a  unit  which,  over  a  small  time 
increment,  has  a  constant  air  concentration  over  its  length 
which  then  serves  as  a  constant  boundary  condition  for 
transport  into  the  fluid  layers. 

Transport  in  the  fluid  and  tissue  layers  presents  a 
much  simpler  analytical  problem.  The  combined  thickness  of 
the  fluid  blanket  and  epithelium  is  small  compared  to  the 
diameter  of  the  airway.  Thus,  for  purposes  of  mass 
transport  modeling,  we  may  neglect  the  cylindrical  nature  of 
the  problem  and  assume  that  the  ring  of  fluid  and  epithelium 
is  stretched  out  along  a  smooth  plane.  This  will  allow  use 
of  trigonometric  functions  rather  than  the  more  complicated 
Bessel  expressions.  Furthermore,  the  upward  movement  of  the 
g lyco-protein  layer  due  to  ciliary  action  is  very  slow  such 
that  the  overwhelmingly  predominant  influence  on  that 
layer's  concentration  is  the  concentration  gradient  due  to 
radial  diffusion,  and  the  convective  component  can  be 
neglected.  Thus,  the  problem  can  be  perceived  as  shown  in 
figure  5.3.  Since,  over  a  small  time  increment,  air 
concentration  in  a  single  airway  is  constant  and  tissue 


u  n  i  form 
cone  C0 


u  n i form 
cone  C  2 


SIMPLIFIED  RIRWRY  MODEL 
FIGURE  5.2 


concentration  is  constant,  we  have  a  one  dimensional 
diffusion  problem  across  four  infinite  planes  separating  the 
air,  fluid,  and  tissue  regions.  Crank**  presents  an 
expression  for  mass  transport  between  two  infinite  planes 
with  constant  boundary  conditions  Cl  and  C2  at  the  two 
planes,  respectively,  and  initial  concentration  profile 
f  (X)  : 

C  =  Cl  +  (C2  -  Cl  )X/L  +  2/tt  [C2  Cos(n7T)  -  Cl]/n  * 

n=l 

Sin(  n  77-X/L)  exp( -DW  n27T2t/L2)  +  2/L  ^  Sin(n7TX/L)  * 

n=l 

exp(-DW  n2  7T2t/L2)/Lf  (X  ’ )  Sin(n  7TX  '/L)dX '  (5.2) 

A) 

where  L  =  distance  between  boundaries 

and  X  =  distance  from  lower  boundary  (X=0). 

Since  the  water  layer  represents  a  reservoir  of  uniform 
concentration  due  to  turbulent  mixing  by  cilia,  we  can  apply 
equation  5.2  across  the  g lyco-protein  layer  with  Cl  equal  to 
the  constant  air  concentration  times  the  air/liquid 
partition  coefficient  and  C2  equal  to  the  uniform  water 
concentration  during  a  small  time  increment.  The  same 
equation  may  be  used  to  describe  transport  in  the  epithelium 
with  Cl  and  C2  equal  to  the  water  concentration  and  tissue 
concentration  respectively.  The  uniform  water  concentration 
must  then  be  adjusted  from  time  increment  to  time  increment 
(t  =  say,  0.01  sec)  based  on  evaluation  of  the  derivative  of 
equation  5.2  at  its  boundaries  and  calculation  of  total  mass 


fluxing  in  or  out  of  the  known  water  volume  during  the  time 


increment.  It  now  remains  to  define  the  air  concentration 
and  tissue  concentration.  Consider  the  beginning  of  an 
inhalation  exposure  at  which  time  a  frontal  wave  of  inspired 
air  enters  the  trachea.  At  some  time  later,  the  wave  has 
moved  into  the  trachea  such  that  the  concentration  behind  it 
is  CO  (inspired  air  concentration)  and  the  concentration 
ahead  of  it  is  0 ,  as  shown  in  figure  5.4.  Since  the  model 
considers  a  constant  air  concentration  throughout  the  airway 
during  a  time  increment,  at  some  point  the  tracheal  air 
concentration  must  be  re-defined  from  0  to  CO.  A  reasonable 
approximation  would  be  to  establish  that  point  at  the  time 
the  wave  extends  half  way  into  the  trachea,  or  one-half 
tracheal  retention  time  at  the  prevailing  flow.  Next,  the 
model  must  re-define  the  concentration  leaving  the  trachea 
(input  concentration  for  next  generation)  based  on  the 
cummulative  mass  absorbed  through  the  wall  as  determined  by 
means  described  above.  This  definition  should  occur  when 
the  wave  reaches  the  end  of  the  trachea,  or  full  retention 
time.  The  next  generation  will  then  use  this  value  to 
define  its  air  concentration  when  half  retention  is  achieved 
for  that  generation.  Wave  entry  times  must  be  established 
for  each  generation  at  the  beginning  of  exposure  and 
redefined  each  time  full  retention  is  reached  for  the 
respective  generation  based  on  the  triangular  flow  function 
over  the  course  of  inhalation  or  exhalation. 


BRONCH I RL  AIR  CONCENTRATION 

FIGURE  5. 4 


The  concentration  leaving  generation  16  is  averaged 
over  the  course  of  an  inhalation  and  is  used  to  define 
alveolar  concentration  (which  is  the  input  concentration 
during  exhalation).  We  assume  that  the  concentration 
entering  alveolated  space  is  instantaneously  equilibrated  in 
blood/gas  exchange  such  that  alveolar  concentration  is 
established  by  a  simple  distribution  over  the  combined 
volume  of  the  respiratory  zone  and  blood  flow.  An 

O 

appropriate  expression  is  given  by  Fi  serova -Be  rgerova  as 


CALV  *  cin  10/(0  +  F 


where 


C  ALV 

Cin 

P/0 


1  + ( F/0 ) PN  (5.3) 

alveolar  concentration 

avg  concentration  leaving  generation  16 

per  fusion/ venti lation  ratio 
1.25  in  resting  man 


PN  =  blood/air  partition  coefficient  (taken  to 
be  equal  to  the  water/air  partition  coef¬ 
ficient  for  the  inhaled  gas) 


Inhaled  gas  or  vapor  concentration  in  the  blood  should 

be  distributed  back  into  lung  tissue  from  the  heart  in  a 

matter  of  seconds  after  undergoing  alveolation.  We  must 

therefore  consider  its  influence  on  lung  tissue 

concentration.  An  expression  for  tissue  concentration  is 

1 5 

taken  from  Ketty  :  (We  assume  here  that  the  tissue/blood 


partition  coefficient  is  1.0.) 


where  Ci  =  tissue  concentration 
Ga  =  PN  CALV 
k  =  F/V 

V  =  lung  tissue  volume  g 

=  600  cc  ( F  i  serova -Be  rgerovaB ) 

F  =  bronchial  tissue  perfusion  rate 
(estimated  at  0.5  liters/min.) 

By  employing  the  model  above  with  input  values  of 
diffusion  coefficient,  partition  coefficient,  exposure 
concentration,  tidal  volume,  and  respiratory  rate,  we  can 
describe  concentrations  and  magnitudes  and  directions  of 
mass  flux  throughout  the  conductive  zone  over  small  time 
increments.  Cumulative  values  can  be  formed  such  that,  at 
the  end  of  a  complete  inhalation/exhalation  cycle,  fractions 
of  inhaled  mass  can  be  calculated  for  alveolar  uptake, 
bronchial  epithelial  uptake,  discharge  to  GI,  and 
exhalation.  Figure  5.5  summarizes  the  function  of  the  model 
and  lists  the  values  ultimately  produced.  As  shown,  the 
fractional  components  of  distribution  should  add  to  1.0  for 
ideal  results.  The  computed  fractional  values  might  be 
expected  to  change  in  proportion  to  each  other  in  subsequent 
breaths  as  blood  and  tissue  concentrations  change  but  should 
stabilize  fairly  quickly  such  that  concentrations  and  rates 
of  flux  remain  the  same  from  end-exhalation  to  end- 
exhalation  . 

Figure  5.6  lists  the  major  assumptions  of  the  model 
described  and  each  is  summarized  below: 


TRANSPORT  MODEL  OVERVIEW 
FIGURE  5.5 


MODEL  ASSUMPTIONS 


1.  Regular  dichotomy 


2.  Bronchial  absorp.  only 


3.  Uniform  air  cone,  prfle 


4  .  Rir/liq  interface  equ  i  1 


5.  Low  concentration 


6.  Slow  ciliary  transport 


FIGURE  5.6 


1.  Weibel's  regular  dichotomy  model  for  lung 
morphology  is  used  which  assumes  symmetrical  bifurcation  of 
airway  branches  with  each  airway  of  a  particuar  generation 
having  identical  dimensions.  Bifurcation  in  the  human  lung 
is  known  to  be  'irregular*. 

2.  The  model  considers  transport  in  the  conductive 
zone  only  and  does  not  deal  with  nasal-pharyngeal  absorption 
or  epithelial  absorption  in  alveolated  airways  of  the 
transitory  and  respiratory  zones.  Thus,  tidal  volume  and 
exposure  concentration  refer  to  air  entering  the  trachea  at 
body  temperature  and  humidity. 

3.  The  solubility  of  the  inhaled  gas  or  vapor  is 
assumed  to  be  such  that  no  significant  air  concentration 
profile  is  developed  in  the  airways  due  to  absorption  at  the 
wall.  Furthermore,  air  flow  patterns  at  bifurcations  are 
such  that  mass  is  uniformly  distributed  across  the  airway 
cross-section  upon  entry  into  an  airway  generation. 

4.  The  air  and  liquid  concentrations  at  the  air/liquid 
interface  are  in  phase  equilibrium  at  the  airway  wall. 

5.  The  inhaled  gas  or  vapor  is  at  low  concentration 
such  that  the  phase  equilibrium  function  is  linear,  being 
described  by  a  partition  coefficient,  and  molecular 
diffusion  can  be  described  using  Fick's  Law. 

6.  The  convective  motion  of  the  mucus  blanket  due  to 
ciliary  action  is  slow  and  is  not  a  significant  factor  in 
mass  transport  of  gases  and  vapors. 


Chapter  6 


RESULTS 

The  model  described  in  chapter  5  was  implemented  using 
fortran  77.  A  source  listing  of  one  version  of  the  main 
routine  (with  comments)  is  provided  in  Appendix  C. 
Significant  cpu  time  was  logged  on  the  Cyber  205  at  Colorado 
State  University  through  the  Institute  for  Computational 
Studies  there  and  the  services  of  Triangle  Universities 
Computation  Center.  Most  of  the  results,  however,  were 
generated  by  the  MAP  6420  array  processor  by  CSP  Inc, 
Billerica,  Mass.  This  system  was  hosted  by  a  VAX  750  at  the 
National  Institutes  of  Health's  National  Biomedical 
Simulation  Resource  facility  at  Duke  University.  Results 
were  cross-checked  using  the  VAX  7  50  with  a  Digital 
Equipment  Corp.  fortran  compiler. 

Several  trial  runs  were  performed  using  various  time 
increments  for  successive  iterations.  It  was  found  that 
reducing  the  time  increment  from  0.01  second  to  0.001  second 
gave  concentration  grid  values  which  differed  by  less  than 
2%.  Subsequent  runs  were  performed  at  an  iteration  exposure 
time  increment  of  0.01  second  to  conserve  cpu  time. 

The  algorithm  was  constructed  to  continue  generating 
data  until  such  time  as  end-exhalation  conditions  differed 
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by  less  than  5%  of  conditions  at  the  most  previous  end- 
exhalation,  as  discussed  in  chapter  5.  Under  this 
criterion,  only  the  data  generated  for  the  last  breath  are 
reported  as  model  results.  in  all  cases,  conditions 
stabilized  within  15  breaths  with  most  stabilizing  in  the 
range  of  3  to  8  breaths,  lower  solubilities  tending  to 
require  longer  stabilization  times. 

Reasonable  ranges  of  solubilities  (PN)  and 
dif fusivities  { DFW )  were  selected  based  on  review  of 
handbook  data  for  hundreds  of  compounds.  However,  with 
higher  solubilities,  model  results  for  some  cases  revealed 
100%  absorption  in  the  conductive  zone.  Since  this  was 
considered  inconsistent  with  a  'relatively  insoluble' 
assumption,  the  solubility  range  was  reduced  and  the  ranges 
studied  were  as  follows: 

[  6.26  x  10-6  <  DFW  <  4.74  x  10-!j  ]  cm2/sec 

[  0.62  <  PN  <  13.0  ] 

where  the  partition  coefficient  (PN)  is  defined  as  described 

in  chapter  4.  Above  a  PN  value  of  13  the  uniform  air 

concentration  assumption  would  not  be  expected  to  hold  with 
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high  diffusivity  values  (  >  4  x  10  cm  /sec).  Under  these 
conditions,  high  bronchial  absorption  would  likely  produce  a 
significant  parabolic  air  concentration  profile.  The 
resulting  lower  air  concentration  at  the  bronchial  walls 
would  limit  further  absorption  and,  thus,  bronchial 


absorption  would  be  lower  and  alveolar  uptake  would  be 

higher  relative  to  model  predictions.  However,  diffusivity 
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values  greater  than  3  x  10  are  rare  and  in  the  range  of  1- 
2  x  10  ^ ,  which  holds  for  most  compounds,  the  PN  range  for 
the  model  may  be  extended  to  approximately  20. 

Operating  within  the  above  ranges  109  cases  were 
studied  with  the  following  critical  input  data: 

Exposure  Concentration  =  375  mg/1 

Tidal  Volume  =  600  cc 

Respiratory  Rate  =  15  breaths/min 

-4 

G-P  layer  thickness  =  7.0  x  10  cm 

-4 

Water  layer  thickness  =  7.0  x  10  cm. 

(Note:  Model  output  is  in  terms  of  concentrations  based  on 

exposure  concentration.  Fractional  partitioning  is  then 
calculated  from  model  output.  Test  runs  at  10,  75,  125,  and 
250  mg/1  exposure  concentration  gives  identical  fractional 
uptake  data,  showing  that  fractional  partitioning  is 
independent  of  concentration.  Thus,  as  one  might 
intuitively  assume,  total  uptake  is  directly  proportional  to 
exposure  concentration.) 

Results  were  analyzed  using  SAS  utilities  (Statistical 
Analysis  System)  to  sort,  plot,  and  compare  data,  and  to  fit 
simple  analytical  expressions  to  some  curves.  Initial  model 
results  from  above  data  are  listed  (sorted  by  PN  and  DFW)  in 
Appendix  D.  The  last  column  in  the  listing  represents  the 
sum  of  the  fractional  compar tmenta 1  uptake  values  which 


should  ideally  add  to  1.0. 

Figure  6.1  is  a  plot  of  model  results  showing  total 
mass  taken  up  by  the  lung  (mass  inhaled  minus  mass  exhaled) 
as  a  fraction  of  mass  inhaled  vs  solubility  (PN).  The  three 
curves  represent  three  values  of  the  diffusion  coefficient, 
the  lower  curve  at  the  bottom  of  the  DFW  range,  the  middle 
curve  at  mid  range,  end  the  top  curve  at  top  of  the  DFW 
range.  The  curves  show  that  solubility  is  the  primary 
influencing  factor  in  total  uptake  while  diffusivity  appears 
to  govern  slightly  the  asymptotic  value  approached  at  higher 
solubilities.  Note  also  that  uptake  is  significant,  even  at 
very  low  solubilities  and  climbs  sharply  up  to  a  PN  value  of 
4,  nearing  an  asymptotic  maximum  beyond  a  PN  of  8. 

Figure  6.2  shows  a  plot  of  the  fraction  taken  up 
through  the  epithelial  walls  of  the  conductive  zone.  Here 
we  see  a  much  greater  dependence  on  diffusivity  as  this  is 
the  primary  transport  pathway  in  which  diffusivity  is  a 
factor.  Diffusivity  may  indeed  play  as  important  a  role 
here  as  solubility.  This  is  illustrated  by  a  least  squares 
non-linear  regression  fit  to  the  plotted  data.  A  reasonable 
fit  is  obtained  with  an  expression  giving  equal  weight  to 
both  variables: 

FREPI  =  1  -  exp(  -3961  DFW  PN )  .  (6.1) 

Figure  6.3  shows  the  same  model  results  for  two  of  the 
diffusivity  values  with  curves  overlayed,  representing 
corresponding  values  predicted  by  equation  6.1.  Appendix  E 


lists  the  complete  model  results  for  fractional  epithelial 
uptake  with  corresponding  values  and  residuals  from  equation 

6.1. 

Figure  6.4  shows  a  plot  of  the  fraction  taken  up  by 
blood/gas  exchange  in  the  respiratory  zone,  again,  for  the 
same  three  diffusivity  values.  Here  we  can  see  evidence  of 
two  terms  which  are  apparently  affecting  the  shape  of  the 
curve.  This  can  easily  be  explained  when  we  consider  that 
alveolar  uptake  is  dependent  not  only  on  solubility  in 
blood/gas  partitioning  but  also  on  the  mass  allowed  to  reach 
the  alveolar  region  from  the  conductive  zone.  This  alveolar 
uptake  rises  with  solubility  up  to  point  where  uptake  is 
limited  by  decreased  input  due  to  rising  absorption  in  the 
conductive  zone.  Regression  analysis  provides  the  following 
empirical  equation  for  alveolar  uptake: 

FRALV  =  exp( -6  28  2  DFW  PN )  -  1/(1  +  348.5  DFW°'5PN)  (6.2) 

Again,  figure  6.5  and  Appendix  F  compare  model  results  with 
values  predicted  by  equation  6.2  as  in  the  case  of  bronchial 
absorption . 

The  fraction  discharged  to  the  GI  tract  by  muco-ciliary 
clearance,  as  shown  in  figure  6.6,  is  not  a  significant 
component  of  distribution.  Unlike  particulate  matter  which 
is  captured  and  accumulated  by  the  mucus  blanket  and  carried 
upward  to  the  trachea,  gas  and  vapor  concentrations  in  the 
fluid  layers  conform  to  the  concentrations  in  adjacent 
regions  without  a  cumulative  build  up.  Thus,  the  fluid 
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layer  concentration  combined  with  a  small  volume  of  fluid 
discharged  in  the  course  of  a  single  exposure  limits 
discharge  to  the  GI  tract  to  extremely  small  values. 

However,  prediction  of  the  discharge  may  be  valuable  in  some 
cases  when  the  impact  of  toxic  substances  on  the  GI  tract  is 
a  concern.  Figure  6.6  shows  that  the  fraction  reaching  the 
GI  tract  is  almost  solely  dependent  on  solubility  and  can  be 
practically  estimated  by  a  linear  regression  according  to 
equation  6.3: 

FRG I  =  5.75  X  10“8PN  (6.3) 

Sensitivity  to  Assumed  Parameters 

The  model,  as  described,  assumes  constant  values  for 
certain  key  physiological  parameters  for  simplicity.  To 
evaluate  the  model's  sensitivity  to  changes  in  these 
parameters,  several  model  runs  were  performed  at  a  constant, 
mid-range  diffusion  coefficient.  For  each  parameter 
considered,  thirteen  runs  were  performed  at  thirteen 
representative  solubility  points  within  the  model  solubility 
range.  Each  run  held  the  parameter  in  question  constant  at 
a  value  substantially  changed  from  the  model  base-line  value 
while  holding  all  other  parameters  constant  at  the  base-line 
value.  Table  6.1  lists  factors  by  which  fractional 
bronchial  epithelial  absorption  and  fractional  alveolar 
uptake  were  changed  from  the  model  base-line  results  at  the 
same  solubility  and  diffusivity  point.  The  factors  listed 
are  averages  of  the  thirteen  values  produced  for  that  case. 


but  none  of  the  calculated  factors  differed  from  their 
average  by  more  than  5%.  The  results  listed  in  table  6.1 
are  discussed  below. 

Respiratory  Rate: 

Increasing  the  respiratory  rate  from  15  to  25  breaths 
per  minute  substantially  reduces  the  fraction  absorbed  in 
the  upper  airways  (0.61  of  base-line  value)  due  to  the 
higher  flow  rate  and  decreased  retention  time  in  this 
region.  Less  absorption  increases  the  mass  entering  the 
alveolar  region;  however,  this  effect  is  offset  somewhat  by 
the  higher  alveolar  ventilation  rate  which  tends  to  lower 
alveolar  uptake.  Hence,  only  a  slight  increase  in 
fractional  alveolar  uptake  is  effected.  Thus,  the  combined 
effect  of  lower  bronchial  absorption  and  only  slightly 
altered  alveolar  uptake  would  appear  to  lower  total 
fractional  uptake.  At  first,  this  appears  to  be  incongruous 
with  increased  respiration.  However,  these  figures  are  for 
uptake  as  a  fraction  of  that  inhaled.  Sinced  a  higher 
respiratory  rate  increases  the  volume  inhaled  per  minute, 
the  total  mass  taken  up  per  minute  would  not  be  expected  to 
decrease  and  may  in  fact  increase. 

Tidal  Volume: 

Increasing  the  tidal  volume  is  precisely  analogous  to 
increasing  the  respiratory  rate  (higher  flow  rate,  smaller 
retention  time).  Thus,  the  same  effect  is  observed.  The 
smaller  factor  for  tidal  volume  relative  to  that  for 


Table  6.1 


Factors  by  which  Baseline  Model  Results  Change 
with  Indicated  Parameter  Change 


PARAMETER  CHANGE 

Increase  resp  rate  to  25 

Double  tidal  volume 

Double  mucus  layer  thickness 

Double  perfusion/volume  ratio 
in  bronchial  tissue 

Half  alveolar  perfusion/ 

ventilation  ratio  (Q/V) 


FRACT.  BRONCH.  FRACT.  ALV 


ABSORP. 

UPTAKE 

0.61 

1.10 

0.47 

1.15 

0.97 

0.98 

0.94 

1.01 

0.80 

0.90 

respiratory  rate  reflects  the  degree  to  which  tidal  volume 
was  altered  compared  to  that  for  repiratory  rate. 

Mucus  Layer  Thickness: 

An  increase  in  the  thickness  of  the  gly co-protein  layer 

may  be  encountered,  particularly  in  the  case  of  upper 

respiratory  disease  or  an  irritant  effect  by  the  exposure 
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gas.  Wanner  points  out  that  the  walls  and  cilia  of  the 
upper  airways  can  support  a  substantially  greater  mucus 
thickness  than  is  normally  found  in  healthy  humans.  For 
this  test,  the  gly co-prote in  thickness  in  the  trachea  was 
doubled  and  linearly  tapered  down  from  generation  to 
generation  to  the  base-line  value  in  generation  16.  Results 
show  only  minor  sensitivity. 

Bronchial  tissue  perfusion/volume  ratio: 

The  perfusion/ volume  ratio  in  bronchial  tissue  is  used 
to  calculate  bronchial  tissue  concentration.  Doubling  this 
ratio  (halving  tissue  volume  or  doubling  perfusion  rate)  has 
only  a  minor  influence  on  model  results.  The  higher  ratio 
produces  a  higher  tissue  concentration  which  slightly  limits 
bronchial  absorption  and  allows  slightly  more  gas  or  vapor 
to  reach  the  alveolar  region. 

Alveolar  perfusion/ventilation  ratio: 

This  test  measured  the  effect  of  reducing  the 
perfusion/ventilation  ratio  in  the  calculation  of  the 
alveolar  equilibrium  air  concentration  only,  apart  from  any 
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change  in  respiratory  rate  or  tidal  volume  (i.e.  a  lowering 
of  the  perfusion  rate).  The  resulting  higher  concentration 
in  arterial  blood  produces  higher  bronchial  tissue 
concentrations  (assuming  the  bronchial  perfusion  rate 
remains  unchanged).  Thus,  bronchial  absorption  is  reduced. 
Alveolar  uptake  also  is  reduced  due  to  the  smaller  blood 
flow . 

Conclusions 

In  summary,  the  present  work  sought  to  develop  a 
simplified  model  of  gas  and  vapor  transport  in  the  human 
lung  to  partition  inhaled  mass  into  compartments  of  alveolar 
blood/gas  exchange,  bronchial  epithelial  absorption,  muco¬ 
ciliary  clearance  to  the  GI  tract,  and  exhalation.  The 
following  statements  summarize  the  results  of  the  work: 

1.  Sufficient  tools  are  available  to  construct  a 
practical  model  of  mass  transport  in  the  lung  and  to  collect 
sufficient  data  for  analysis. 

2.  Total  uptake  by  the  lung  is  significant  even  at 
very  low  solubilities.  (fig  6.1) 

3.  Uptake  by  direct  blood/gas  exchange  in  the 
respiratory  zone  decreases  with  solubility  after  peaking  at 
a  relatively  low  solubility  due  to  bronchial  absorption  in 
the  upper  airways.  (fig  6.4) 

4.  Diffusivity  in  the  fluid  layers  and  epithelium  of 


the  bronchial  wall  is  a  significant  factor  in  bronchial 
absorption,  perhaps  as  much  as  solubility.  (fig  6.3) 


5.  Muco-c il iary  clearance  of  inhaled  gas  or  vapor  to 
the  GI  tract  is  insignificant  as  a  fraction  of  mass  inhaled. 

6.  Sufficient  model  results  display  a  clear  trend  in 
uptake  patterns  allowing  empirical  development  of  equations 
estimating  fractional  compartmental  uptake  as  a  function  of 
solubility  and  diffusivity  (eq'ns  6.1  thru  6.3).  These 
equations  may  then  be  employed  in  expanded  modeling  efforts. 

7.  In  the  absence  of  experimental  data,  sufficient 
methods  exist  for  the  estimation  of  solubilities  and 

d if f usivities  using  readily  available  handbook  data  for  most 
gases  and  vapors  (chapter  4). 
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APPENDIX  A 


LAMINAR  FLOW  VELOCITY  PROFILE 
IN  A  CIRCULAR  TUBE 


Expressions  describing  the  velocity  distribution  in 

flow  through  various  geometries  are  described  in  Bird, 

2 

Stewart,  and  Lightfoot  .  For  the  case  of  laminar  flow  in  a 
circular  tube,  consider  a  cylinder  of  length  L  and  radius  R 
as  shown  in  figure  A.l.  We  begin  by  assuming  that  ideal 
laminar  flow  exists,  i.e.  flow  everywhere  is  parallel  to  the 
tube  wall  and  entrance  and  exit  perturbations  in  flow 
streamlines  are  absent.  We  consider  an  elemental  volume 
within  the  tube  as  a  hollow  shell  of  thickness  Ar  and 
length  L  and  seek  to  write  an  equation  expressing  a  momentum 
balance  over  this  volume.  Since  a  parabolic  velocity 
distribution  is  developed  due  to  frictional  resistance  at 
the  wall,  a  'shearing  force'  per  longitudinal  area  may  be 
thought  of  along  the  surface  of  the  shell  which  we  will 
denote  by  '7".  The  rate  of  momentum  transfer  (force)  across 
the  shell  may  then  be  expressed  in  terms  of  this  shearing 
force  per  area: 


Rate  of  momentum  in  =  2  IT  rL  7"|  (A.l) 

Rate  of  momentum  out  =  2  TT  rL  7“|  „ .  .  (A. 2) 

'  r+  A  r 


In  the  z  direction,  the  rate  of  momentum  (force  or  mass- 
velocity  per  time)  can  be  expressed  in  terms  of  the  mass 
entering  or  leaving  and  its  velocity: 

Rate  of  momentum  in  =  2  IT  r  A  rv  (  f*  v )  I  _  (A. 3) 

'  Z  =  (J 

Rate  of  momentum  out  -  2tt  r  6rv(  j5  v)  I  _  (A. 4) 

z  =  L 

The  sum  of  these  rates  of  momentum  transfer  over  the 
elemental  volume  must  equal  the  total  force  acting  on  the 
element.  These  forces  can  be  described  as: 


Gravity  force  = 

27Tr  ArLfg 

(A. 5) 

Pressure  force  = 

2  v  r  ArP(0  ) 

(A. 6) 

at  z=0 

Pressure  force  = 

2  7T  r  ArP  (  L) 

(A. 7) 

at  z=L 

The  sum  of  equations  A. 1  thru  A. 7  must  be  zero: 

2  77"  rL  7^|  r  -  2  7TrL  7-|  r+  Ar  +  2  7T  r  A  r  P  v2  |  z=Q  - 

2 

2  7T  r  A  r  Z5  v  |  r  +  2  77’r  4rLfq  + 

'  z=L 

27T  r  A  r  (P  (0  )  -  P  (  L )  1  =  0  (A. 8) 

We  now  simplify  the  equation  by  assuming  the  compression  of 
the  air  over  the  length  of  the  tube,  given  small  static 
pressure  differences  within  the  system,  is  negligible.  Thus 


87 


we  may  consider  the  velocity  at  z=0  to  equal  the  velocity  at 
z=L  and  the  third  and  fourth  terms  above  cancel.  If  we  then 
divide  by  2  TT  L  A  r  and  take  the  limit  as  Ar  approaches  zero, 
we  have: 


d/dr(  r  7“*)  =  [P  ( 0  )  -  P(L)]r/L  +  Pgr  (A. 9) 

Recognizing  that  gravity  is  acting  in  the  same  manner  at  z=0 
and  z=L,  we  may  drop  the  right  most  term  by  assuming  that 
P(0)  and  P(L)  represent  the  combined  effect  of  static 
pressure  and  gravity.  Upon  integration  we  then  have 

7"  =  LMO)  -  P  (  L)  ]  r/2L ,  (A. 10) 

the  constant  of  integration  being  zero  to  satisfy  conditions 
at  r=0. 

The  shearing  force  per  area  may  also  be  described  by 
Newton's  law  of  viscosity  as 

'7~  =  -/dv/dr  (A. 11) 

Combining  equations  A. 10  and  A. 11,  we  have 

dv/dr  =  ~[P(0)  -  P(L)]r/2yuL 


and  integrating, 
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v  =  -  tP  (0  )  -  P  ( L)  ]  rz/4//L  +  C  (A. 12) 

Sice  the  velocity  is  zero  at  r=R 

C  =  [P(0)  -  P  ( L)  ]  R2/4yUL  (A. 13) 


and 


V  =  [P  ( 0  )  -  P  ( L)  ]  R2  [1  -  (  r/R)  2]  /4/if  L  (A. 14) 

To  simplify  this  equation,  one  may  use  the  above  expression 
for  v  and  calculate  the  average  velocity  over  the  cross- 
section  by  integrating  [vrdrd#  ]  from  0  to  R  and  from  0  to 
2  and  dividing  by  the  cross-sectional  area.  This  procedure 
will  show  that  the  average  velocity  is 

VAVG  =  tP(0)  ~  P(L)]  r2/8//L  (A. 15) 

Thus  the  velocity  profile  may  be  written  in  terms  of  the 
average  velocity  as 


v  =  2  v..7„[l  -  (r/R)2] 


(A. 16) 


APPENDIX  B 


EQUATIONS  OF  CHANGE  FOR 
TURBULENT  AIR  FLOW 

Referring  to  figure  3.1,  we  desire  to  write  a  mass 
balance  over  the  elemental  volume  for  the  case  of  turbulent 
flow.  Here  we  consider  the  mass  flux  in  the  z  direction  to 
consist  of  a  convective  flux,  which  is  an  average  flux  over 
time,  plus  a  turbulent  flux,  which  represents  variations 
from  the  time  averaged  value  due  to  turbulence. 

Nz  =  vC  +  v'C'  (B.l) 

where  v  is  the  time  averaged  velocity  at  r.  Likewise,  since 
the  turbulent  flux  is  expected  to  act  radially  as  well,  the 
flux  in  the  r  direction  consists  of  a  diffusive  flux  plus 
the  turbulent  flux. 

Nr  =  -Da  be/  dr  +  v'C’  (B.2) 

For  conditions  in  the  turbulent  core  of  the  stream,  we  use 
the  von  Karman  Similarity  Hypothesis  to  describe  the 
turbulent  flux  as  follows:  (see  Bird,  Stewart,  and 
Lightfoot^ ) 


Nz  =  vC  -  K2j(  d  v/  d  r)2 /(  d2v/  ^r2)2|  dc/  d  z 


( B .  3  ) 


Nr  =  -da  dc/  6  r  -  K2  |  (  dv/  dr)2/ (  d2v/  dr2)2  \  dc/  dr 

( B .  4  ) 

where  K  =  0.36  or  0.40  depending  on  the  investigation 
referenced . 

Experimental  investigations  have  shown  the  velocity 
profile  in  a  circular  tube  for  turbulent  flow  to  be  roughly 

v  =  1.25  vAVG[l  -  r/RG]0,143  (B .  5 ) 

Taking  the  first  and  second  derivative  of  the  above 
expression  with  respect  to  r,  we  may  then  substitute  in 
equations  B.3  and  b.4,  giving  values  of  Nz  and  Nr  which  may 
now  be  used  to  evaluate  equation  3.3  for  the  turbulent  case 

-  d  Nz/  dz  -  <)(Nr  r)/r  dr  =  dc /  d  t  (3.3) 

The  resulting  equation  for  mass  transport  in  the  turbulent 
core  is 

Da[  d2C/  dr2  +  d  C/r  dr]  +  vftVG  ( 1  -  r/RG)0,143  * 

[0.032  RG  (1  -  r/RG)  ^2C/  d z2  -  1.25  dc/  d z  =  0.032  RG 
(1  -  r/RG)  d2 C/  ^r2  +  0.032  RG/r  *  (1  -  r/RG)  dc/  dr  - 


-0.037  dc/  dr]  =  dc/  d  t 


(B.6) 


Near  the  wall  of  the  tube  (outside  the  turbulent  core), 

motion  is  dampened  by  resistance  at  the  wall,  and  the 

turbulent  flux  can  be  described  by  Deissler's  Empirical 
9 

Formula  so  that  the  fluxes  are  written  as 

Nz  =  vC  -  n2vr[l  -  exp(  -n2vr/ )  ]  d  C/  dz  (B.7) 

Nr  =  -D  dc/  dr  -  n2vr[l  -  exp(  -n2vr/\T  )]  d  C/  dr  (B.8) 

where  V*  =  kinematic  viscosity 

and  n  =  0.124  (empirical  constant). 

Using  equation  B.5  for  velocity  and  substituting  equations 
B.7  and  B.8  into  equation  3.3,  one  arrives  at  the  mass 
transport  equation  for  the  near-wall  region.  The  resulting 
long  expression  is  not  reproduced  here. 

The  location  of  the  boundary  separating  the  turbulent 
core  from  the  near-wall  region  is  a  function  of  the  medium's 
velocity  and  viscosity  and  can  be  determined  in  terms  of  the 
dimensionless  value  S+: 

S+  =  0.36  S2{  d  v/  d  S)//  ( B .  9  ) 


Using  equation  B.5  for  velocity  and  recognizing  that 
r  =  RG  -  S,  the  above  equation  can  be  re-written  as 


S+  =  0.064  vAVG(RG  -  r  r  *14  3/(  V'RG0  *143) 


(B.10) 


Values  of  r  for  which  S+  is  less  than  26  are  considered  to 
be  in  the  near-wall  region. 


APPENDIX  C 

FORTRAN  SOURCE  LISTING 


The  following  pages  display  a  FORTRAN-77  source  listing 
for  the  transport  model  described  in  chapter  5.  Only  the 
main  routine  showing  the  general  model  algorithm  is  listed. 
More  detailed  calculations  from  sub-routines  are  not  shown. 


C  <*  Model  lung  for  insoluble  case  *) 

PROGRAM  INSOL 

C**********************************************************) 
C*  *) 

C*  This  program  performs  a  running  calculation  of  the  *) 

C*  concentration  of  a  particular  species  (vapor  or  gas)  *) 
C*  in  various  regions  of  the  human  lung.  It  is  intended  *) 
C*  to  approximate  the  case  of  a  relatively  insoluble  *) 

C*  vapor  or  gas  in  water  such  that  no  significant  concert-  *) 
C*  tration  profile  is  developed  within  a  single  *) 

C*  bronchial  segment.  It  relies  on  constant  input  values  *) 
C*  of  physiological  parameters  and  species  *) 

C*  characterist ics.  *) 

C*  Weibel’s  morphological  model  of  the  human  lung  is  used.*) 
C*  The  inner  surface  of  each  bronchiole  is  assumed  to  be  *) 
C*  bathed  with  a  smooth  continuous  bi layer  fluid  of  *) 

C*  viscous  glycoprotein  on  tc<D  and  non-vi scous  water  *) 

C*  beneath.  A  characterist ic  thickness  is  assumed  for  *) 
C*  the  epithelia  separating  the  fluid  from  well  perfused  *) 
C*  tissue.  The  transport  equation  employed  is  from  Crank  *) 
C*  (Eqn  A. 16),  'The  Mathematics  of  Diffusion’.  *) 


C*  PROGRAMMER: 

C*  FUNCTIONS  CALLED: 
C* 


C*  PROCEDURES  CALLED: 
C* 


M.  L.  SHELLEY 


RAVG 

RINTGR 

FLOW 

RETNTN 

TRMIII 

TERMIV 

DERVTV 


PRTCON 

CHKCND 


C**********************************************************) 


CHARACTERS®  SPEC 

(*  exposure  gas  or  vapor  *) 

REAL  PNCOEF 

(*  water/air  partition  coef  *) 

REAL*8  DFCOFW 

(*  diffusion  coef  of  species  *) 
(*  in  water  *) 

REAL  ATCOEF 

(*  activity  coef  for  solubility  in  water  *) 

REAL  PS 

(*  saturation  vapor  pressure  *) 


REAL  VC 


c 


REPL*8  EXPCON 


(*  critical  volume  *) 


vv 

.  * .  ■» 

kV.> 


■§$ 


t  - 


tv 


C 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


REPL  MPXTIM 


<*  exposure  cone  (micro-grams/ec)  *) 


c 

<* 

REPL*8 

TRDIP 

c 

REPL*8 

TRLEN 

c 

REPL*8 

PTRDIP 

c 

REPL*8 

PTRLEN 

c 

REPL*8 

TIDVOL 

c 

REPL*8 

GLYC0 

c 

<* 

REPL*8 

GLYC16 

c 

<* 

REPL*8 

WPTR0 

c 

<* 

REPL*8 

WPTR16 

c 

<* 

REPL*B 

EPIT0 

c 

<* 

REPL*8 

EPIT16 

c 

<* 

(*  dimensions  of  trachea  *) 

(*  "  *) 

<*  ernperical  value  for  gen  4  to  16  *) 
( *  "  * ) 

(*  tidal  volume  of  lung  *) 


*> 

*) 

*) 

*) 

*) 


REPL  DELTIM 

<*  time  increment  for  iteration  *> 

REPL  PCTSTP 

(*  acceptable  signif  for  lung  stability  *> 
REPL*8  GLYVEL 

(*  mucus  transport  rate  *> 

INTEGER  MPXPNT 

<*  number  of  points  defining  cone  prfle  *) 
INTEGER  RESRPT 

(*  respiratory  rate  *) 

INTEGER  MXTERM 

<*  terms  for  series  approx irnat ions  *) 
CHPRPCTER*8lZl  PSTRSK 
REPL*8  PI 

REPL*8  PI RCJN  <0: 16, £) 

<*  bronchial  air  concent rat  ion  *) 

REPL  ENTIME  <0:16) 

<*  time  from  beginning  of  exposure  at  *) 
<*  which  the  frontal  wave  enters  a  *> 
(*  bronchiole  *> 

REPL*8  PMPSTR (0:16) 

<*  cum  mass  transported  from  air  to  *) 
(*  glyco  throughout  a  generation  *) 

<*  during  one  retention  time  *) 

REPL*8  FLUXTO<0: 16, 4) 

<*  current  flux  across  interface  *) 
REPL*8  NFLXTO (0: 16, 4) 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


<*  flux  across  interface  at  next  time  incr  *) 
REAL*8  MASSTO (0: 16, 4) 

(*  mass  transferee!  across  interface  in  *) 
<*  an  individual  bronchiole  in  one  #) 

<*  time  increment  *) 

REAL *8  LNGCON (0 : 16, £, 30) 

<*  cone  prfle  in  lung  region  *) 

REAL*8  LNGC01 (17) 


(*  for  MAP  application  only  *) 


INTEGER  JMAP,  KMAP 


(*  "  *) 
REAL*8  NLNCON<0: 16, £, 30) 

<*  cone  one  time  increment  later  *) 
REAL*8  WATCON (0:16) 


REAL*8 

REAL*8 

REAL*8 

REAL*8 

REAL*8 


(*  cone  in  water  layer  *) 


NWTCON (0:16) 

(*  cone  one  time  increment  later  *) 


LNGDIM(0: 16, 6) 

(*  dimensions  of  lung  *) 

Cl (0:50) 

<*  avg  cone  discharged  to  resp  in  breath  *) 
UPTAKE <0:50) 

<*  avg  uptake  rate  over  breath  (mass/sec)  *) 
GIDSCR (0:50) 


(*  avg  discharge  rate  to  GI  over  breath,  m/s  *) 
REAL*8  EPDSCR (0:50) 

(*  avg  discharge  rate  thru  epi  wall  *) 


INTEGER  GLY,  WAT,  EPI,  TIS,  GLYCOP,  EPITHL,  DIA,  VOL, 
1  LEN,  GLYTH, 

1  WATTH,  EPITH,  TOGI,  TOEXH,  TOTIS,  STRD,  IN, 

1  OUT 

<*  array  subscripts  *) 

REAL*8  TLCGLY(30),  TLCEPI (30) ,  NTCGLY (30) 

<*  temporary  arrays  for  parameter  passing  *) 


REAL  CURT I M 

C  <*  elapsed  exposure  time  *) 

REAL*8  ALVCON 

C  (*  alveolar  concentration  *) 

REAL*8  SYSCON 

C  <*  systemic  concent rat  ion  *) 

REALMS  MASTOG 

C  <*  mass  transfered  to  GI  during  current  time  increment  *) 
REAL*8  CMSTOG 

C  <*  cumrnulative  mass  discharged  to  GI  *) 

REALMS  CMSTOR 

C  <#  cum  mass  passed  to  respiratory  zone  *) 

C  (*  in  current  breath  *) 

REAL *8  MASTOE 

C  <*  mass  passed  to  exhalation  in  current  time  incr  *) 


REAL*8  CMSTOE 

C  <*  cum  mass  exhaled  in  current  breath  *) 


nn  n  n  n  n  n  no 


REAL*8  CMSTOT 


<*  cum  mass  transfered  to  tissue  *) 

REAL  TMBRTH 

(*  time  from  beginning  of  exposure  at  which  *) 
<*  the  current  inhalation  or  exhalation  began  *) 
REAL  CUMVOL 

<*  cummulative  volume  for  calc  initial  entry  times  *) 
INTEGER  I 

<*  loop  counter  *) 

INTEGER  GEN 

<*  generation  number  (trachea  =  0)  *) 

INTEGER  POINT 

(*  evaluation  point  for  cone  function  *) 
INTEGER  ITERTN 

(■*  iteration  counter  *) 

INTEGER  BREATH 

(*  number  of  full  breath  (inh  +  exh)  from  *) 
(*  beginning  of  exposure  *) 

LOGICAL  INHLTN 

(*  states  whether  inh  or  exh  *) 

LOGICAL  STBLZD 

<*  states  condition  of  luno  *) 


REAL *8  EPS I 


<*  zero  variable  *) 


C**************************^******************************) 

C  INITIALIZE  CONSTANTS 

C ********************************************************* ) 


DATA 

SPEC/’ 111’ 

/ 

DATA 

ATCOEF, 
/0.  1, 

PS,  EXPCON,  MAXTIM, 

0.01,  3.  75D+02,  60.05, 

TRDIA 

1. 8D+00/ 

DATA 

TRLEN 
/ 12. 0D+00/ 

DATA 

PTRDIA, 
/l. 3D+00, 

PTRLEN,  TIDVOL,  GLYC0, 

2.  5D+00,  6. 0D+02,  7. 0D-04, 

GLYC16 

7. 0D-04/ 

DATA 

WATR0 
/7. 0D-04/ 

DATA 

WATR16, 
/7. 0D-04, 

EPIT0,  EPIT16,  DELTIM 

3.  0D-03,  3. 0D-03,  0.01, 

,  PCTSTA 
0.  05/ 

DATA 

GLYVEL 
/8. 3D-03/ 

DATA 

MAXPNT, 

RESRAT,  MXTERM,  VC 

1  / 15,  15,  12,  40.0/ 


C**********************************************************) 


DATA  PI/3. 1415926536D+00/ 

DATA  EPSI/1. 0D-06/ 

DATA  IN,  OUT,  GLY,  WAT,  EPI,  TIS,  GLYCOP,  EPITHL,  DIA 


1  /l,  2,  1, 

2,  3,  4, 

1, 

-  - - f 

DATA  VOL,  LEN 

1  /£,  3/ 

DATA  GLYTH,  WATTH, 

EPITH,  TOGI, 

TOEXH, 

TOT IS, 

STRD 

1  /  4,  5, 

fc,  1, 

£, 

3, 

4/ 

C*********************************************************) 
C*  <*  MAIN  PROGRAM  *>  *) 
C*********************************************************) 


C  (*  Initialize  program  variables  *> 

PNCOEF  =  1243/ (ATCOEF  *  PS) 

DFCOFW  *  2. 27D-04/ <0. £85  *  VC**1 . 048) **0.  6 
SYSCON  =  0 
DO  110  GEN  =  0, 16 
AIRCON(GEN, IN)  =  0 
AIRCON (GEN, OUT)  =  0 
AMASTR (GEN)  =  0 
WATCON (GEN)  =  0 
DO  120  POINT  =  1 , MAX PNT 

LNGCON ( GEN, GLYCOP, POINT)  =  0 
LNGCON<GEN, EPITHL, POINT )  =  0 
120  CONTINUE 

FLUXTO(GEN, GLY)  =  0 
FLUXTO(GEN, WAT)  =  0 
FLUXTO (GEN, EPI )  =  0 
FLUXTO (GEN, TIS)  =  0 
110  CONTINUE 


C 


70 


(*  Calculate  rnorpholog ica  1  dimensions  of  lung  *) 


DO  70  GEN  =  0,  3 
LNGDIM (GEN,  DIA) 
LNGDIM (GEN, LEN) 
LNGDIM (GEN, VOL) 

1 


TRDIA  *  EXP<-0. 388  *  BEN) 
TRLEN  *  EXP (-0.920  *  GEN) 
PI/4  *  (LNGDIM (GEN, DIA) > **£ 
LNGDIM (GEN, LEN)  *  2. 0**6EN 


* 


CONTINUE 

DO  80  GEN  =4,16 
LNGDIM (GEN, DIA) 
1  (0.0062  *  GEN  - 

LNGDIM (GEN, LEN) 
LNGDIM (GEN, VOL) 


=  PTRDIA  *  EXP (GEN  * 

0.  293)  ) 

=  PTRLEN  *  EXP (-0.17  *  GEN) 

=  PI/4  *  (LNGDIM (GEN, DIA) ) **£  * 
LNGDIM (GEN, LEN)  *  2. 0**GEN 


1 


99 


80  CONTINUE 


C 


90 


(*  Calculate  layer  thicknesses  *) 


DO  90  GEN  =0,16 
LNGDIM (GEN, GLYTH) 

LNGDIM (GEN, WATTH) 

LNGDIM (GEN, EPITH) 

CONTINUE 


GLYC0  -  GEN  * 

(GLYC0  -  GLYC16) /16 
WATR0  -  GEN  * 

(WATR0  -  WATR16) /16 
EPIT0  -  GEN  * 

(EPIT0  -  EPIT16J/16 


C  (*  Initialize  all  entry  times  *) 

ENT I  ME (0)  =  0 
CUMVOL  =  0 
DO  100  GEN  =  1,16 

CUMVOL  =  CUMVOL  +  LNGDIM (GEN  -  l,VOL> 

ENTIME (GEN)  =  SORT (450  *  CUMVOL/ (RESRAT**2  * 
1  TIDVOL)) 

100  CONTINUE 


Cl (0)  =  0. 0 
UPTAKE  (0)  =0.0 
EPDSCR  (0)  =0.0 
GIDSCR (0)  =  0.  0 


ALVCON 

= 

0 

CURTIM 

= 

0 

INHLTN 

= 

.  TRUE. 

CMSTOG 

= 

0 

CMSTOR 

= 

0 

CMSTOE 

= 

0 

CMSTOT 

0 

TMBRTH 

= 

0 

BREATH 

= 

1 

STBLZD 

= 

.  FALSE. 

WRITE  (ASTRSK,  ’  <c'A40>’> 

’  **************************************** 

’  **************************************** 

MASTOE 

= 

0 

C  (*  Write  program  heading  *) 

WRITE (6, 510)  ASTRSK,  SPEC,  ASTRSK,  ASTRSK 

WRITE (6, 511)  EXPCON,  DFCOFW,  PNCOEF, 

1  MAXPNT,  MAXTIM,  TIDVOL,  RESRAT,  DELTIM, 

1  MXTERM, 

1  PCTSTA,  TRDIA,  TRLEN,  ATCOEF,  PS,  VC 

510  FORMAT (///IX, A80/1X, 8X, ’ RESULTS  OF  LUNG  TRANSPORT  ’, 

1  ’ MODEL  FOR  ’  , 

1  ’INSOLUBLE  SPECIES  ’, 

1  A20/1X, A80///1X, A80//1X,  38X,  ’  CASE’  /  //> 

511  FORMAT ( 1 X, 8X, 

1  ’EXP  CONC  =  ’  ,  D8.  2,  6X,  ’  DIFF  COEF  =’ ,  D9.  £,  IX, 

1  ’  PARTITION  COEF  =  ’ , 


M7,£X, 


1  E8.  £//lX, £X, ’ PROFILE  POINTS  * 

1  ' TOTAL  EXP  TIME  =  ’ , F7. £, 

1  4X,  ’  TIDAL  VOLUME  =  ’  ,  D8.  £//’  RESPIRATORY 

1  » ROTF  s  »  T7 

1  ’TIME  INCREMENT’*  »,F8.  4,  £X, 

1  ’  SERIES  TERMS  =  ’,17// 

1  ’  STABILITY  SIGNIF  =  ’ , F7. £, 5X, 

1  ’TR  DIAMETER  =  ’,D7.2,7X, 

1  ’TR  LENGTH  =  ’ , D7. £// 

1  ’  ACTIVITY  COEF  =  ’ , F7. £, IX, 

1  ’SAT  VAPOR  PRESS  =  ’ ,  F7.  £,  IX, 

1  ’CRITICAL  VOLUME  =  ’,F7.£/> 


ITERTN  =  0 

(*  Begin  exposure  and  cummulative  calculation  *> 

>  CONTINUE 

ITERTN  =  ITERTN  +  1 

(*  Begin  single  iteration  *) 

DO  £90  GEN  =  0,16 

(*  If  time  to  re-assign  concentration  in  *) 

IF  (CURTIM  .GE.  < ENT I ME (GEN)  + 

1  RETNTN <’ HALF’ ,LNGDIM (GEN, VOL), 

1  TMBRTH, ENTIME (GEN) , TRDI A, TRLEN, PTRDI A 

1  PTRLEN, RESRAT, TIDVOL) ) ) 

1  THEN 

IF  ( INHLTN) 

1  THEN 

IF  (GEN  .NE.  0) 

1  THEN 

AIRCON (GEN,  IN)  *  AIRCON (GEN  -  1,0UT) 
ELSE 

AIRCON (GEN, IN)  =  EXPCON 
END  IF 
ELSE 

<*  NOT  INHALATION  *) 

IF  (GEN  .NE.  16) 

1  THEN 

AIRCON (GEN, IN)  =  AI RCON (GEN  +  l,OUT> 
ELSE 

AIRCON (GEN, IN)  =  ALVCON 
END  IF 

END  IF 

<*  if  inhalation  *) 

END  IF 

(*  if  time  to  re-aasign  *) 


(*  If  time  to  re-assign  cone  out  and  entry  time  *) 
IF  (CURTIM  .GE.  (ENTIME (GEN)  + 

1  RETNTN (’ FULL’,  LNGDIM (GEN,  VOL), 


1 

1 

1 


THEN 


TMBRTH, ENTIME (GEN) , TRDIA, TRLEN, PTRDI A, 
PTRLEN, RESRAT, TIDVOL) ) ) 


C  (*  Calculate  concentration  out  *) 

AIRCON  (GEI1  HUT)  =  < AI RCON (GEN, IN)  * 

1  LNGDIM (GEN, VOL)  - 

1  AMASTR (GEN) ) /LNGDIM (GEN,  VOL) 

IF  <AIRCON(GEN,  OUT)  .  LT.  0.0) 

1  THEN 

AI RCON (GEN, OUT)  =  0. 0D+00 
END  IF 

C  (*  Re-init  ial  ize  cum  mass  transferee!  *) 

AMASTR (GEN)  =  0 

C  (*  Re-assign  entry  time  *) 

ENTIME (GEN)  =  CURTIM 

END  IF 


C  (*  Create  temporary  array  variable  for  pararn  passing  *) 

DO  310  I  =  1 , MAXPNT 

TLCGLY ( I )  ■  LNGCON (GEN, GLYCOP, I ) 

TLCEPI(I)  =  LNGCON (GEN, EPITHL,  I) 

310  CONTINUE 

C  (*  Calculate  next  f(r)  points  for  glyco  and  epi  regions  *) 
DO  300  POINT  =  1, MAXPNT 

IF  ( (AIRCQN(GEN, IN)  . LT.  EPSI)  .AND.  (WATCON (GEN) 

1  .LT.  EPSI)) 

1  THEN 

NLNCON(GEN, GLYCOP,  POINT)  =  0 
ELSE 

NLNCON (GEN, GLYCOP, POINT)  =  PNCOEF  * 

1  AI RCON (GEN, IN)  + 

1  (WATCON (GEN)  -  PNCOEF  *  AI RCON (GEN, IN) )  * 

1  POINT/ (MAXPNT  +  1)  + 

1  TRMIII (POINT,  AIRCON(GEN, IN)  * 

1  PNCOEF,  WATCON (GEN)  ,  LNGDIM (GEN,  GLYTH)  , 

1  MAXPNT, 

1  MXTERM, DFCOFW, DELTIM,  PI)  + 

1  TERMIV (POINT, TLCGLY, LNGDIM (GEN,  GLYTH)  , 

1  MAXPNT, MXTERM, DFCOFW,  DELTIM,  PI) 

END  IF 

IF  ( (WATCON (GEN)  . LT.  EPSI)  .AND.  (SYSCON  . LT. 

1  EPSI)) 

1  THEN 

NLNCON (GEN, EPITHL, POINT)  =  0 
ELSE 

NLNCON (GEN, EPITHL, POINT)  =  WATCON (GEN)  + 
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END  IF 


(SYSCON  -  WATCON (GEN) )  *  POINT/ 
(MAXPNT  +  1)  + 

TRMI 1 1 ( POINT, WATCON (GEN) ,  SYSCON, 
LNGDIM (GEN, EPITH) ,  MAXPNT, 

MXTERM, DFCOFW, DELTIM, PI)  + 

TERMI V ( POINT, TLCEPI , 

LNGDIM (GEN, EPITH) ,  MAXPNT, 

MXTERM, DFCOFW, DELTIM,  PI > 


CONTINUE 


<*  Calculate  fluxes 
IF  ( <AIRCQN(GEN,  IN) 
.LT.  EPSI )  ) 
THEN 

NFLXTO (GEN, GLY) 
NFLXTO (GEN, WAT) 
ELSE 

NFLXTO (GEN, GLY) 


between  regions  *) 

.LT.  EPSI)  .AND.  (WATCON (GEN) 


NFLXTO (GEN, GLY)  *  -DFCOFW  *  DERVTV  (’ LOWER’ , 
TLCGLY, AIRCON (GEN, IN)  *  PNCOEF, 

WATCON (GEN) , LNGDIM (GEN, GLYTH) , 

MAXPNT, MXTERM, DFCOFW, 

DELTIM, PI) 

NFLXTO (GEN, WAT)  *  -DFCOFW  *  DERVTV  <’ UPPER’ , 
TLCGLY, AIRCON (GEN,  IN)  *  PNCOEF, 

WATCON (GEN) , LNGDIM (GEN, GLYTH) , 

MAXPNT, MXTERM, DFCOFW, 

DELTIM, PI) 

END  IF 


IF  ( (WATCON (GEN)  .LT.  EPSI)  .AND.  (SYSCON  . LT. 

EPSI) ) 

THEN 

NFLXTO (GEN, EPI)  =  0 

NFLXTO(GEN, TIS)  =  0 
ELSE 

NFLXTO (GEN, EPI)  =  -DFCOFW  *  DERVTV (’ LOWER’ , 
TLCEPI, WATCON (GEN) ,  SYSCON, 

LNGDIM (GEN, EPITH) ,  MAXPNT, 

MXTERM, DFCOFW, DELTIM, PI ) 

NFLXTO (GEN, TIS)  =  -DFCOFW  *  DERVTV (’ UPPER’ , 
TLCEPI, WATCON (GEN) , SYSCON, 

LNGDIM (GEN, EPITH) , MAXPNT, 

MXTERM, DFCOFW, DELTIM, PI > 

END  IF 

<*  Calculate  mass  transferee!  between  regions  *) 
MASSTO ( GEN, 6L Y )  ■  (NFLXTO (GEN, GLY)  + 

FLUXTO (GEN, GLY) ) /£  * 

DELTIM  *  PI  *  LNGDIM (GEN, LEN)  *  LNGDIM (GEN, DIA) 
MASSTO (GEN, WAT)  =  (NFLXTO (GEN, WAT)  + 


V'- 


n  n 


103 


1  FLUXTO (GEN,  WAT) > /£  * 

1  DELTIM  *  PI  *  LNGDIM (GEN, LEN)  *  LNGDIM (GEN, DIP) 

MASSTO (GEN,  EPI )  =  (NFLXTO (GEN, EPI )  + 

1  FLUXTO (GEN,  EPI) ) /£  * 

1  DELTIM  *  PI  *  LNGDIM (GEN, LEN)  *  LNGDIM (GEN, DIA) 

MASSTO (GEN, TIS)  »  (NFLXTO (GEN,  TIS)  + 

1  FLUXTO (GEN, TIS)) /£  * 

1  DELTIM  *  PI  *  LNGDIM (GEN,  LEN)  *  LNGDIM (GEN,  DIA) 
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C 


CONTINUE 


320 


(*  end  GEN  loop  *) 


(*  Update  concentration  and  curnmulat  i  ve  values  *> 

(*  Create  ternp  array  variable  for  param  passing  *) 
DO  320  I  =  1 , MAXPNT 

TLCGLY ( I )  =  LNGCON(0, GLYCOP, I) 

NTCGLY  < I )  =  NLNCON(0, GLYCOP, I) 

CONTINUE 

MASTOG  =  ( RAVG (TLCGLY, MAXPNT, AI RCON (0, IN)  *  PNCOEF, 
WATCON (0) )  +  RAVG (NTCGLY, MAXPNT, 

AIRCON(0, IN)  *  PNCOEF, (WATCON (0)  *  PI  * 
LNGDIM (0,  DIA)  *  LNGDIM (0,  LEN)  * 


LNGDIM (0, WATTH) 
MASSTO (0,  WAT)  - 
LNGDIM (0,  DIA)  * 
LNGDIM (0,  LEN)  * 
TRDIA  * 

LNGDIM <0,GLYTH) 


MASSTO (0,EPI) )/ (PI  * 

LNGDIM (0, WATTH) )>)/2  *  PI  * 
*  GLYVEL  *  DELTIM 


CMSTOG  =  CMSTOG  +  MASTOG 


IF  ( INHLTN) 

1  THEN 

CMSTOR  = 

1 

ELSE 

MASTOE  = 

1 

1 

CMSTOE  = 
END  IF 


CMSTOR  +  AI RCON (16, OUT)  *  DELTIM  * 
FLOW (CURTIM, TMBRTH, RESRAT, TIDVOL) 

AI RCON (0, OUT)  *  DELTIM  * 

FLOW (CURTIM, TMBRTH, 

RESRAT, TIDVOL) 

CMSTOE  +  MASTOE 


DO  330  GEN  =0,16 


CMSTOT  =  CMSTOT  +  MASSTO (GEN, TIS)  *  2. 0**GEN 


AMASTR (GEN)  =  AMASTR (GEN)  +  MASSTO (GEN, GLY)  * 
1  2. 0**GEN 


WATCON (GEN)  =  (WATCON (GEN)  *  PI  *  LNGDIM (GEN, DIA)  * 

1  LNGDIM (GEN, LEN)  *  LNGDIM (GEN, WATTH)  + 

1  MASSTO (GEN, WAT)  -  MASSTO (GEN, EPI ))/( PI 
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*  LNGDIM (GEN,  DI A)  *  LNGDIM (GEN, LEN)  * 
LNGDIM (GEN, WATTH) ) 

IF  (WATCON (GEN)  .  LT.  0.0) 

THEN 

WATCON (GEN)  =  0. 0D+00 
END  IF 

DO  340  POINT  =  1 , MAXPNT 

LNGCON (GEN, GLYCOP, POINT)  = 

NLNCON (GEN, GLYCOP,  POINT) 

LNGCON (GEN, EPITHL, POINT)  = 

NLNCON (GEN, EPITHL,  POINT) 

CONTINUE 

FLUXTO(GEN, GLY)  =  NFLXTO (GEN, GLY) 

FLUXTO(GEN, WAT)  =  NFLXTO (GEN, WAT) 

FLUXTO (GEN, EPI )  =  NFLXTO (GEN, EPI ) 

FLUXTO (GEN, TIS)  =  NFLXTO (GEN, TIS) 


CONTINUE 


<*  end  GEN  Iood  *) 


<*  Increment  time  *) 

CURTIM  =  CURTIM  +  DELTIM 

<*  If  end  of  inhalation  or  exhalation  *) 

IF  (CURTIM  .GE.  (TMBRTH  +  30/RESRAT) ) 

1  THEN 

TMBRTH  =  CURTIM 

DO  370  GEN  =  0, 16 

ENT I ME (GEN)  =  CURTIM 
CONTINUE 

IF  ( INHLTN) 

1  THEN 

ALVCON  =  CMSTOR/TIDVOL  *  (1/(1  +  l.£5  * 

1  PNCOEF) ) 

SYSCON  =  PNCOEF  *  ALVCON  *  (1  -  EXP (-0.0139  * 

1  CURTIM)) 

Cl (BREATH)  =  CMSTOR/TIDVOL 
CMSTOR  =  0.0 

INHLTN  =  .FALSE. 

CALL  PRTCON ( ITERTN,  CURTIM,  ASTRSK, 

1  BREATH,  INHLTN,  TMBRTH, 

1  SYSCON,  LNGCON,  WATCON,  AIRCON,  GLYCOP, 

1  EPITHL, 

1  IN,  ALVCON,  CMSTOR,  CMSTOE,  CMSTOT,  CMSTOG) 

ELSE 

<*  exhalation  *) 
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UPTAKE (BREATH)  =  (EXPCON  *  TIDVOL  -CMSTOE)  * 
1  RESRAT/60 

CMSTOE  =0.0 

GIDSCR (BREATH)  =  CMSTQG  *  RESRAT/60 
CMSTOG  =  0.0 

EPDSCR (BREATH)  =  CMSTOT  *  RESRAT/60 
CMSTOT  =  0.0 

INHLTN  =  .TRUE. 

CALL  CHKCND  <CI, UPTAKE, STBLZD, EPDSCR, GIDSCR, 

1  ASTRSK, BREATH, CURTIM, PCTSTA, EPSI ) 

BREATH  =  BREATH  +  1 

END  IF 

C  (*  if  inhalation  else  exhalation  *) 

END  IF 

C  (*  If  end  inhalation  or  exhalation  *) 


C  (*  end  iteration  loop  *> 

IF  ((CURTIM  .LT.  MAXTIM)  .AND.  (.NOT.  STBLZD))  THEN 
GO  TO  £70 
END  IF 

C  (*  Write  program  heading  *) 

WRITE (6, 512)  ASTRSK,  SPEC,  ASTRSK,  ASTRSK 

WRITE<6, 513)  EXPCON,  DFCOFW,  PNCOEF, 

1  MAXPNT,  MAXTIM,  TIDVOL,  RESRAT,  DELTIM, 

1  MXTERM, 

1  PCTSTA,  TRDIA,  TRLEN,  ATCOEF,  PS,  VC 

512  FORMAT (///IX, A80/ IX, 8X, ’ RESULTS  OF  LUNG  TRANSPORT  ’, 

1  ’MODEL  FOR  INSOLUBLE  SPECIES  ’, 

1  A20/1X, A80///1X, A80//1X,38X, ’CASE’ ///) 

513  FORMAT (IX, 8X, 

1  ’EXP  CONC  =  ’  ,  D8.  2,  6X,  ’  DIFF  COEF  =’,D9.£,  IX, 

1  ’PARTITION  COEF  =  ’, 

1  E8.  2//lX,£X,  ’  PROFILE  POINTS  =  ’,I7,2X, 

1  ’TOTAL  EXP  TIME  =  ’ , F7. 2, 

1  4X,* TIDAL  VOLUME  =  ’ , D8. £// 

1  ’  RESPIRATORY  RATE  =  ’,I7,£X, 

1  ’TIME  INCREMENT  =  ’ ,  F8.  4,  £X, 

1  ’  SERIES  TERMS  =  ’ , 17// 

1  ’  STABILITY  SIGNIF  =  ’ ,  F7. 2,  5X, 

1  ’TR  DIAMETER  =  ’,D7.£,7X, 

1  ’TR  LENGTH  =  ’,D7.£// 

1  ’  ACTIVITY  COEF  =  ’,F7.£, IX, 

1  ’SAT  VAPOR  PRESS  =  ’ ,  F7.  2,  IX, 

1  ’CRITICAL  VOLUME  =  ’,F7.£/) 

CALL  PRTCON ( ITERTN,  CURTIM,  ASTRSK,  BREATH,  INHLTN, 

1  TMBRTH, 


SYSCON,  LNGCON,  WATCON,  AIRCON, 
GLYCOP,  EPITHL, 

IN,  ALVCON,  CMSTOR,  CMSTOE,  CMSTOT 
CMSTOG) 

(*  List  stabilization  parameters  *) 

WRITE (6, 389)  ASTRSK 

FORMAT ( IX, A80/1X, ’ BREATH’ , 9X, ’ Cl’ , 9X, ’ UPTAKE  RATE’ , 5X 
’EPI  DSCH  RATE’ , 3X, ’ GI  DSCH  RATE’//) 

DO  380  I  =  0, BREATH-1 

WRITE (6, 390)  I,  Cl < I ) ,  UPTAKE (I),  EPDSCR(I), 
GIDSCR(I) 

FORMAT  <3X,  13,  5X,  D13.  7,  3X,  D13.  7,  3X,  D13.  7,  3X,  D13.  7/) 
CONTINUE 

(*  Write  condition  of  normal  job  termination  *> 

IF  (CURTIM  .GE.  MAXTIM) 

THEN 

WRITE<6, 350)  CURTIM 

FORMAT (/////////IX, ’TIME  =  ’,F8.4, 

’  ******  NORMAL  TERMINATION’/) 

ELSE 

I DUMMY  =  BREATH  -  1 

WRITE <6, 241)  ASTRSK,  IDUMMY,  CURTIM,  Cl ( I DUMMY) , 
UPTAKE < I DUMMY ) ,  EPDSCR ( IDUMMY) , 

GIDSCR( IDUMMY) 

FORMAT(////1X,A80//1X,30X, ’STABILITY  PARAMETERS’ / 1 X, 
31 X,  ’END  OF  ’,  14, 

’  BREATHS’ /IX, 28X, 

’EXPOSURE  TIME  =  ’ , F8. 4///1 X, 24X, 

’  CNC  DISCHARGING  TO  ALV  *  ’ , D1 1 . 4// 1 X, 16X, 

’  UPTAKE  PER  SECOND  = 

’, Dll. 4/ /IX, 20X, 

’  DISCHARGE  RATE  TO  TISSUE  =  ’ , D1 1 . 4// 1 X, 26X, 
’DISCHARGE  RATE  TO  GI  ■  ’,D11.4//) 

WRITE (6, 360)  CURTIM,  BREATH 

FORMAT (/////////IX, ’CONDITIONS  STABLE  AT  ’,F8.4, 

’  (BREATH  ’, 14,’ )’, 

’  ******  NORMAL  TERMINATION’/) 

END  IF 

STOP 

END 


(*  MAIN  PROGRAM  *) 


APPENDIX  D 


TABULATED  MODEL  RESULTS 


The  following  pages  list  model  results  for  109 
hypothetical  gases  or  vapors  (OBS) .  The  columns  are  defined 
as  follows: 


PN  =  Partition  coefficient  (as  defined  in  chapter  4) 

£ 

DFW  =  Diffusion  coefficient  in  water  (cm  /sec) 

BRTH  =  Number  of  breaths  required  to  reach  stable 
condit ions 

FRUPTK  =  Fraction  of  inhaled  mass  taken  up  <1  -  fract 

exhaled) 

FREPI  =  Fraction  of  inhaled  mass  taken  up  through 

epithelial  absorption  in  the  conductive  zone 

FRALV  =  Fraction  of  inhaled  mass  taken  up  by  alveolar 
blood/gas  exchange 

FRGI  =  Fraction  of  inhaled  mass  taken  up  by  muco-ciliary 
discharge  to  the  GI  tract 

CHK  =  (1  -  FRUPTK)  +  FREPI  +  FRALV  +  FRGI 
(ideally  equal  to  1.0) 


BS  PN 


DFW  BRTH  FRUPTK  FREPI  FRftLV  FRBI 


CHK 


1 

0.617 

6. 26E-6 

2 

0.617 

1 . 03E-5 

3 

0.617 

1.72E-5 

4 

0.617 

2. 85E-5 

5 

0.617 

4. 74E-5 

6 

0.  885 

7.  16E-6 

7 

0.885 

8. 49E-6 

8 

0.  885 

1.07E-5 

9 

0.  885 

1. 57E-5 

10 

1.080 

7.  16E-6 

11 

1.080 

8.49E-6 

12 

1.080 

1.07E-5 

13 

1.080 

1.57E-5 

14 

1.  100 

7. 16E-6 

15 

1.  100 

8. 49E-6 

16 

1.  100 

1.07E-5 

17 

1.  100 

1. 57E-5 

18 

1. 310 

6.  26E-6 

19 

1. 310 

1. 72E-5 

20 

1.310 

4. 74E-5 

21 

1. 340 

7. 16E-6 

22 

1.340 

8. 49E-6 

23 

1.  340 

1.07E-5 

£4 

1.340 

1 . 57E-5 

£5 

1.400 

7. 16E-6 

26 

1.400 

8. 49E-6 

27 

1.400 

1.07E-5 

£8 

1. 400 

1.57E-5 

29 

1. 490 

7. 16E-6 

30 

1.490 

8. 49E-6 

31 

1.490 

1.07E-5 

32 

1. 490 

1.57E-5 

33 

1.560 

6. 26E-6 

34 

1.560 

1 . 72E-5 

35 

1.56 

4. 74E-5 

36 

1.75 

7. 16E-6 

37 

1.  75 

8. 49E-6 

38 

1.  75 

1. 07E-5 

39 

1. 75 

1.57E-5 

40 

1.  81 

7. 16E-6 

41 

1.81 

8.49E-6 

42 

1. 81 

1. 07E-5 

43 

1.81 

1. 57E-5 

44 

2.  01 

7. 16E-6 

45 

2.01 

8.49E-6 

46 

£.01 

1.07E-5 

47 

£.01 

1.57E-5 

48 

£.  £5 

7.  16E-6 

49 

£.  £5 

8. 49E-6 

50 

£.  £5 

1.07E-5 

si 

£.  £5 

1. 57E-5 

52 

2.35 

7.  16E-6 

5 

0.54 
0.53 
0.44 
5 
5 


1.01E-7 
1.03E-7 
1.03E-7 
1 . 03E-7 
1 • 

I. 39E-7 
1.39E-7 
1.39E-7 
ii.39E-7 
1 . 62E-7 
1 . 64E-7 
1 . 64E-7 
1.64E-7 
I. 65E-7 
1.67E-7 
1 . 67E-7 

1 . 96E-7 
1 . 96E-7 
1.78E-7 
96E-7 
96E-7 
96E-7 
96E-7 
96E-7 
95E-7 
96E-7 
96E-7 
13E-7 
13E-7 
13E-7 
13E-7 
S. 13E-7 
2. 13E-7 
2.  13E=7 

iQC-7 


1.00 

1.00 

1.01 

1.00 

0.  99 
0.99 
0.  99 
1. 01 
0.99 
0.99 
0.99 
1. 01 
0.  98 
0.  98 
1.01 
1. 01 
0.99 
0.  99 
0.99 
9 
9 


0.97 

1.00 

0.98 


0.  59  3.  £0E— 7  0.  98 


3.  3 

3._3' 
3.43 
3.  43 
3.  43 
3.  43 
3.43 
3.57 
3.57 
3.  57 
3.57 


4.  11 
4.  11 
4.  1 1 
4.  1  1 
4.  1 1 


5.  lu 
5.  1£ 
5.  1£ 


74 

7.  16E-6 

74 

8.  49E-6 

74 

1. 07E-5 

1. 7£E-5 
4.  74E-5 


6.  £6E-6 
1.  03E-5 
1. 7£E-5 
£. 85E-5 
4. 74E-5 
7. 16E-6 
8.  49E-6 
1. 07E-5 
1 


6. £6E-6 
1. 03E-5 
1. 7£E-5 
£. 85E-5 
4. 47E-5 


7. 16E-6 
8.  49E-6 
1.07E-5 
1. 57E~C 


95 

5.  8£ 

6.  £6E-6 

96 

5.  8£ 

1. 7£E-5 

97 

5.  8£ 

4.  47E-5 

98 

6.  99 

6. £6E-6 

99 

6.99 

1. 7£E-5 

100 

6.99 

4. 74E-5 

101 

8.  30 

6.  £6E-6 

1 0£ 

8.  30 

1 . 7£E-5 

103 

8.  3 

4.  74E-5 

104 

9.  79 

6.  £6E-6 

0.  59 

4. 80E-7 

0.  95 

0.  50 

4.  80E-7 

0.97 

0.  £9 

4. £7E-7 

0.  96 

0.  59 

4. 98E-7 

0.  95 

0.  56 

4. 98E-7 

0.  96 

0.  50 

4. 98E-7 

0.  97 

0.  40 

4. 6£E-7 

0.  97 

0.  £7 

4. 45E-7 

0.  96 

0.  57 

6. ££E-7 

0.  94 

0.  54 

6.  ££E-7 

0.94 

0.  53 

6.  05E-7 

0.  95 

0.47 

6.  05E-7 

0.96 

0.  56 

6.  93E-7 

0.93 

0.  43 

6. 58E-7 

0.  95 

0.  18 

6.  ££E-7 

0.93 

0.  53 

8.  17E-7 

0.91 

0.  39 

7.  8£E-7 

0.  94 

0.  13. 

7.  £9E-7 

0.91 

0.  50 

9.  60E-7 

0.  90 

0.  35 

9. £5E-7 

0.  93 

0.09 

8.71E-7 

0.  89 

0.  47 

1.  l£E-6 

0.  88 
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OBS 

PN 

DFW 

BRTH 

FRUPTK 

FREPI 

FRALV 

FRGI 

CHK 

105 

9.  79 

1.72E-5 

3 

0.88 

0.  51 

0.  30 

1. 06E-6 

0.  92 

106 

9.79 

4.  74E-5 

3 

0.94 

0.74 

0.07 

1.01E-6 

0.87 

107 

1£.30 

6. 26E-6 

3 

0.  80 

0.25 

0.  41 

1. 38E-6 

0.86 

108 

12.30 

1. 72E-5 

3 

0.  89 

0.  57 

0.  £4 

1. 33E-6 

0.  92 

109 

12.  30 

4.  74E-5 

3 

0.95 

0.81 

0.  05 

1.26E-6 

0.91 

PFREPI 

0.  015 
0.  025 
.  041 
0.  067 
0.  109 
0.  025 
0.  029 
0.  037 
0.  054 
0.  030 
.  036 
0.  045 
0.  065 
0.  031 
0.  036 
0.  046 
0.  066 
0.  032 
0.  085 
0.  218 
0.  037 
0.  044 
0.  055 
0.  080 
0.  039 
0.  046 
0.  058 
0.  083 
0.  041 
0.  049 
0.  061 
0.  089 
0.  038 
0.  100 
0.  254 
0.  048 
0.  057 
0.  071 
0.  103 
0.  050 
0.  059 
0.  074 
0.  106 
0.  055 
0.  065 
0.  082 
0.  117 
0.  062 
0.  073 
0.  091 
0.131 
0.  064 


PRESID 

-0. 01956 
-0. 01975 
-0. 02200 
-0. 02731 
-0. 03550 
-0. 02228 
-0. 02157 
-0. 02111 
-0. 02064 
-0. 02295 
-0. 02210 
-0. 02102 
-0. 02047 
-0. 02286 
-0. 02193 
-0. 01986 
-0. 01996 
-0. 02438 
-0. 01740 
— 0. 01625 
~0. 02266 
-0. 02169 
-0. 01959 
-0. 01704 
-0.02280 
-0. 02097 
-0. 01819 
-0. 01681 
-0. 02308 
-0. 02115 
-0. 01824 
-0. 01668 
-0. 02457 
-0. 01267 
-0. 00315 
-0. 02219 
-0. 01942 
-0. 01708 
-0. 01 186 
-0. 02223 
-0. 01938 
-0. 01693 
-0. 01 156 
-0. 02088 
-0. 01888 
-0. 01395 
-0. 00617 
-0. 02042 
-0. 01629 
—0. 01 263 
-0. 00410 
-0. 02022 


2.  35 

8. 49E-6 

0.  06 

0.  076 

-0. 01596 

2.  35 

1. 07E-5 

0.  08 

0.  095 

-0.  01210 

2.  35 

1. 57E-5 

0.  13 

0.  136 

-0.  00328 

£ .  51 

7.  16E--6 

0.  05 

0.  069 

-0. 01844 

2.  51 

8. 49E--6 

0.  07 

0.  081 

-0. 01569 

2.  51 

1. 07E-5 

0.  09 

0.  101 

-0.  00863 

£.  51 

1. 57E-5 

0.  14 

0.  145 

-0. 00252 

£.74 

7. 16E-6 

0.  06 

0.  075 

-0. 01783 

£.  74 

8. 49E-6 

0.  08 

0.  088 

-0. 01217 

£.  74 

1. 07E-5 

0.  10 

0.  110 

-0.  00731 

£.  74 

1. 57E-5 

0.  16 

0.  157 

0. 00473 

£.  77 

6. 26E-6 

0.  05 

0.  066 

-0.  0£0£7 

£.  77 

1. 7£E-5 

0.  18 

0.  172 

0.  00700 

£.  77 

6. £SE-6 

0.  45 

0.  406 

0.  0484£ 

3.  3i? 

6. £6E-6 

0.  06 

0.  079 

-0.  01946 

o  ■ 

1 . 7£E-5 

0.  ££ 

0.  £01 

0. 01878 

3.  30 

4. 74E-5 

0.  51 

0.  46£ 

0.  04606 

3.  37 

7. 16E-6 

0.  08 

0.  091 

-0.  01412 

3.  37 

8. 49E-6 

0.  10 

0.  107 

-0.  01027 

3.  37 

1. 07E-5 

0.  13 

0.  133 

-0.  00033 

3.  37 

1. 57E-5 

0.  £0 

0.  189 

0.  00892 

3.  43 

6. 26E-6 

0.  06 

0.  08£ 

-0.  01935 

3.  43 

1. 03E-5 

0.  13 

0.  131 

-0. 00242 

3.  43 

1. 72E-5 

0.  £3 

0.  £08 

0.01933 

3.  43 

£. 85E-5 

0.  36 

0.  321 

0.  04289 

3.  43 

4. 74E-5 

0.  5£ 

0.  475 

0.  04480 

3.  57 

7. 16E-6 

0.  08 

0.  096 

-0.  01375 

3.  57 

8. 49E-6 

0.  1 1 

0.  113 

-0.  00603 

3.  57 

1. 07E-5 

0.  14 

0.  140 

0.  00046 

3.  57 

1. 57E-5 

0.  ££ 

0.  199 

0.  01718 

3.  9£ 

6. 26E-6 

0.  08 

0.  093 

-0.  01645 

3.  9£ 

1. 72E-5 

0.  £6 

0.  £34 

0.  02168 

3.  9£ 

4. 47E-5 

0.  56 

0.  521 

0.  0398£ 

4.  1  i 

6. £6E-6 

0.  08 

0.  097 

-0.  01607 

4.  1 1 

1. 03E-5 

0.  16 

0.  154 

0.  00169 

4.  1 1 

1. 72E-5 

0.  £7 

0.  £44 

0.  02275 

4.  11 

£. 85E-5 

0.  41 

0.  371 

0.  04311 

4.  1 1 

4. 47E-5 

0.  58 

0.  538 

0.03781 

5.  1  £ 

7.  16E-6 

0.  13 

0.  135 

-0. 00795 

5.  1£ 

8. 49E-6 

0.  16 

0.  158 

0.  00241 

5.  1  £ 

1. 07E-5 

0.  £1 

0.  195 

0.  01089 

5.  12 

1. 57E-5 

0.  30 

0.  £73 

0.  03058 

5.  8£ 

6. £6E— 6 

0.  13 

0.  134 

-0.  00703 

5.  B£ 

1. 72E-5 

0.  36 

0.  3£7 

0. 03343 

5.  8£ 

4. 47E-5 

0.  67 

0.  665 

0.  00558 

6.  99 

6.  £6E— 6 

0.  15 

0.  159 

-0.00711 

6.  99 

1. 72E-5 

0.  41 

0.  379 

0.  03103 

6.  99 

4. 74E-5 

0.  71 

0.  731 

-0.  02285 

8.  30 

6. £6E-6 

0.  18 

0.  186 

-0.  00235 

8.  30 

1. 72E-5 

0.  46 

0.  43£ 

0. 02625 

8.  3 

4. 74E-5 

0.  73 

0.  789 

-0. 05781 

9.  79 

6. £6E-6 

0.  £1 

0.  £16 

-0.  00467 

OBS 

PN 

DFW 

FREPI 

PFREPI 

PRESID 

105 

9.  79 

1. 72E-5 

0.  51 

0.  487 

0.01878 

1 06 

9.  79 

4. 74E-5 

0.74 

0.841 

-0. 09294 

107 

12.30 

6. 26E-6 

0.25 

0.  2S3 

-0. 01025 

108 

12.  30 

1. 72E-5 

0.57 

0.  567 

0. 00293 

109 

12.30 

4. 74E-5 

0.  81 

0.901 

-0. 09547 

APPENDIX  F 


ALVEOLAR  UPTAKE 
(MODEL  RESULTS  vs  PREDICTED) 


The  following  pages  list  model  results  for  109 
hypothetical  gases  or  vapors  (OBS)  with  corresponding 
predicted  values  from  an  empirically  derived  expression. 


PN 

DFW 

FRALV 

PFRALV 


Partition  coefficient  (as  defined  in  chapter  4) 

2 

Diffusion  coefficient  in  water  (cm  /sec) 

Fraction  of  inhaled  gas  or  vapor  taken  up  by 
alveolar  blood/gas  exchange 


Predicted  FRALV 
exp (-6282  DFW  PN) 


0.  5 

1/(1  +  348.5  DFW  PN) 


PRESID  =  Predicted  residual 


116 


OBS 

PN 

DFW 

FRflLV 

PFRALV 

PRESID 

1 

0.  617 

6.  £6E-6 

0.  400 

0.  386 

0. 07418 

£ 

0.617 

1.03E-5 

0.  397 

0.  369 

0. 0£816 

3 

0.617 

1. 7£E-5 

0.  39£ 

0.407 

-0.01489 

4 

0.617 

£. 85E-5 

0.  376 

0.430 

-0. 05384 

5 

0.617 

4. 74E-5 

0.  376 

0.  4£9 

-0. 05£99 

6 

0.  885 

7.  16E-6 

0.  456 

0.  413 

0. 04£89 

7 

0.  885 

8. 49E-6 

0.  453 

0.  4£7 

0. 0£6 1 3 

8 

0.885 

1. 07E-5 

0.  451 

0.  444 

0. 006££ 

9 

0.  885 

1.57E-5 

0.  467 

0.466 

0. 000£9 

10 

1. 080 

7.  16E-6 

0.  493 

0.  454 

0. 03898 

1 1 

1. 080 

8. 49E-6 

0.  491 

0.  467 

0. 0£359 

1£ 

1. 080 

1.07E-5 

0.  485 

0.  48£ 

0. 00356 

13 

1.080 

1.57E-5 

0.  501 

0.  498 

0. 00378 

14 

1.  100 

7.  16E-6 

0.  491 

0.458 

0. 03£59 

15 

1.  100 

8. 49E-6 

0.  488 

0.  471 

0. 01736 

16 

1.  100 

1. 07E-5 

0.  509 

0.  485 

0. 0£4£8 

17 

1.  100 

1.57E-5 

0.499 

0.  500 

-0. 0015£ 

18 

1. 310 

6.  £6E-6 

0.531 

0.  483 

0. 04769 

19 

1. 310 

1. 7£E-5 

0.499 

0.  5££ 

-0. 0£37£ 

£0 

1. 310 

4.  74E-5 

0.  445 

0.436 

0.00971 

£1 

1. 340 

7.  16E-6 

0.  5£5 

0.  497 

0. 0£836 

££ 

1.  340 

8.  49E-6 

0.  5£3 

0.  507 

0.01 5£5 

£3 

1.340 

1.07E-5 

0.517 

0.  518 

-0. 00089 

£4 

1. 340 

1. 57E-5 

0.  5£8 

0.  5£5 

0. 00£64 

£5 

1.400 

7. 16E-6 

0.  549 

0.  505 

0.04411 

£6 

1.  400 

8. 49E-6 

0.544 

0.515 

0. 0£890 

£7 

1.400 

1. 07E-5 

0.  539 

0.  5£5 

0. 01370 

£8 

1. 400 

1.57E-5 

0.  5£3 

0.  530 

-0. 00743 

£9 

1. 490 

7. 16E-6 

0.  544 

0.  517 

0. 0£734 

30 

1. 490 

8. 49E-6 

0.  539 

0.  5£6 

0. 01300 

31 

1.490 

1.07E-5 

0.  531 

0.  534 

-0. 00345 

3£ 

1. 490 

1. 57E-5 

0.515 

0.536 

-0. 0£159 

33 

1.560 

6. £6E-6 

0.  541 

0.517 

0. 0£453 

34 

1. 560 

1. 7£E-5 

0.  5£5 

0.  538 

—0. 0 1 ££9 

35 

1.56 

4. 74E-5 

0.  440 

0.  418 

0. 0££4£ 

36 

1.75 

7. 16E-6 

0.  563 

0.  544 

0. 01 83£ 

37 

1.75 

8.49E-6 

0.  560 

0.551 

0. 009££ 

38 

1. 75 

1. 07E-5 

0.  555 

0.555 

-0. 00045 

39 

1. 75 

1.57E-5 

0.539 

0.  549 

-0. 01010 

40 

1.  81 

7.  16E-6 

0.  568 

0.  550 

0. 018£3 

41 

1. 81 

8. 49E-6 

0.  563 

0.556 

0. 00707 

4£ 

1.81 

1.  07E-5 

0.  555 

0.559 

-0. 00433 

43 

1.81 

1. 57E-5 

0.  541 

0.  551 

-0. 00940 

44 

£.01 

7. 16E-6 

0.  579 

0.  566 

0. 01 30£ 

45 

£.01 

8. 49E-6 

0.  571 

0.  570 

0. 00117 

46 

£.01 

1. 07E-5 

0.  563 

0.570 

—0.  007 1  £ 

47 

£.  01 

1. 57E-5 

0.544 

0.555 

-0. 01130 

48 

£.  £5 

7.  16E-6 

0.  584 

0.  581 

0. 00303 

49 

£.  £5 

8. 49E-6 

0.  579 

0.  58£ 

-0. 00380 

50 

£.  £5 

1.07E-5 

0.  571 

0.  579 

-0. 00845 

51 

£.  £5 

1.57E-5 

0.  549 

0.557 

-0.00815 

5£ 

£.  35 

7.  16E-6 

0.587 

0.  586 

0. 0003£ 

DFW 

FRALV 

PFRALV 

PRES  ID 

8. 49E-6 

0.  581 

0.587 

-0. 00555 

1. 07E-5 

0.  57 1 

0.582 

-0.01139 

1.57E-5 

0.549 

0.  558 

-0. 00821 

7. 16E-6 

0.  589 

0.  594 

-0. 00455 

8.  49E-6 

0.  584 

0.593 

-0. 00891 

1.07E-5 

0.  571 

0.586 

-0.01509 

1.57E-5 

0.  547 

0.  557 

-0. 01011 

7. 16E-6 

0.592 

0.  603 

-0. 01075 

8. 49E-6 

0.  587 

0.  600 

-0. 01297 

1.07E-5 

0.  573 

0.  589 

-0.01593 

1. 57E-5 

0.  544 

0.554 

-0. 01013 

6. 26E-6 

0.  597 

0.  604 

-0. 00665 

1.72E-5 

0.  533 

0.  541 

-0.00813 

6. 26E-6 

0.  360 

0.  303 

0. 05250 

6. 26E-6 

0.  600 

0.  620 

-0. 02038 

1. 72E-5 

0.520 

0.527 

-0. 00673 

4. 74E-5 

0.  325 

0.  262 

0. 06317 

7.  16E-6 

0.  592 

0.  618 

-0. 02594 

8. 49E-6 

0.  581 

0.  609 

-0. 02800 

1.07E-5 

0.565 

0.  591 

-0. 02541 

1. 57E-5 

0.  533 

0.  540 

-0. 00693 

6. 26E-6 

0.  600 

0.  623 

-0. 02323 

1. 03E-5 

0.  571 

0.  594 

-0. 02352 

1. 72E-5 

0.517 

0.522 

-0. 00510 

2. 85E-5 

0.432 

0.406 

0. 02637 

4. 74E-5 

0.  317 

0.  252 

0. 06559 

7. 16E-6 

0.  592 

0.621 

-0. 02865 

8. 49E-6 

0.  581 

0.  610 

-0. 02907 

1. 07E-5 

0.  563 

0.589 

-0. 02672 

1. 57E-5 

0.  523 

0.  535 

-0.01188 

6. 26E-6 

0.  595 

0.631 

-0.03611 

1. 72E-5 

0.  504 

0.  505 

-0. 00067 

4. 47E-5 

0.285 

0.  215 

0. 07024 

6. 26E-6 

0.  592 

0.  633 

-0. 04058 

1. 03E-5 

0.  557 

0.588 

-0. 03047 

1 . 72E-5 

0.  496 

0.497 

-0. 00130 

2. 85E-5 

0.  400 

0.  363 

0. 03658 

4. 47E-5 

0.272 

0.  202 

0. 06999 

7. 16E-6 

0.  565 

0.621 

-0. 05578 

8. 49E-6 

0.  544 

0.  600 

-0.05571 

1. 07E-5 

0.  525 

0.  563 

-0. 03720 

1. 57E-5 

0.  472 

0.  480 

-0. 00758 

6.  2SE-6 

0.  560 

0.631 

-0. 07080 

1. 72E-5 

0.  432 

0.  427 

0. 00507 

4. 47E-5 

0.  178 

0.  110 

0. 06823 

6.  26E-6 

0.  533 

0.619 

-0. 08536 

1.  72E-5 

0.  389 

0.  380 

0. 00954 

4. 74E-5 

0.  129 

0.068 

0. 06087 

6. 26E-6 

0.  501 

0.  600 

-0. 09876 

1. 72E-5 

0.  347 

0.  331 

0.01577 

4. 74E-5 

0.  090 

0.  037 

0. 05350 

6. 26E-6 

0.  467 

0.576 

-0. 10891 

OBS 


PN 


DFW 


FRALV  PFROLV 


PRES ID 


105 

9.  79 

1 . 72E-5 

0.  298 

0.261 

0.01722 

106 

9.79 

4. 74E-5 

0.  067 

0.013 

0. 05332 

107 

1£.  30 

6. 26E-6 

0.  411 

0.531 

-0. 12025 

106 

12.  30 

1.72E-5 

0.235 

0.211 

0.02401 

109 

12.  30 

4. 74E-5 

0.050 

0.  000 

0. 05725 

v**  k '/«  ■.**-  ;  > -% 

k  «-  viiAlAj Julii  'j 


I'  ^  Vn*  dt  *.» 


APPENDIX  G 


The  following  tables  list  equations  and  coefficients 

used  in  the  prediction  of  activity  coefficients  and  hydrogen 

bonding  factors  for  use  in  determining  the  partition 

coefficient  (PN)  as  described  in  chapter  4.  The  tables  are 

2 1 

adapted  from  Reid,  Prausnitz,  and  Sherwood  and  from 

„  ^  ,13 

Hayduk 


TABLE  G.l 


Correlating  Constants  for  the  Calculation 
of  Activity  Coefficients  in  Water 
(adapted  from  Reid,  Prausnitz,  and  Sherwood^1) 


(All  values  are  for  25  degrees 

C  unless 

otherwise 

stated . ) 

SOLUTE 

J 

K 

M 

EQ'N 

n-Ac  ids 

-1  .000 

0  .622 

0.490 

(a) 

n-Primary  Alcohols 

-0.995 

0.622 

0.558 

(a) 

n-Secondary  Alcohols 

-1.220 

0.622 

0.170 

(b) 

n-Tertiary  Alcohols 

-1.740 

0.622 

0.170 

(c) 

Alcohols,  General 

-0.525 

0.622 

0.475 

(d) 

n-Allyl  Alcohols 

-1.180 

0.622 

0.558 

(a) 

n-Aldehydes 

-0.780 

0.622 

0.3  20 

(a) 

n-Alkene  Aldehydes 

-0 .7  20 

0.622 

0.3  20 

(a) 

n-Ketones 

-1.475 

0.622 

0.500 

(b) 

n-Acetals 

-2.556 

0.622 

0.486 

(e) 

n-Ethers  (20  deg  C) 

-0.770 

0.640 

0.195 

(b) 

n-Nitr iles 

-0.587 

0.622 

0.760 

(a) 

n-Alkene  Nitriles 

-0.5  20 

0.622 

0.760 

(a) 

n-Esters  (20  deg  C) 

-0.930 

0.640 

0.260 

(b) 

n-Formates  (20  deg  C) 

-0.58  5 

0.640 

0.260 

(a) 

n-Monoalkyl  Chlorides 
(20  deg  C) 

1.265 

0.640 

0.073 

(a) 

n-Paraffins  (16  deg  C) 

0.688 

0.642 

0 

(a) 

n-Alkyl  Benzenes 

3.5  54 

0.622 

-0.466 

(f ) 

TABLE  G . 1  (cont. ) 


( a )  LOG  y 

(b)  LOG  y 

( c )  LOG  t 

(d)  LOG  t 

(e)  log  y 

(f )  log  y 

Where 

N'  ,N"  ,N'  ' 


EQUATIONS 


J  +  K  (  N )  +  M/N 

J  +  K(N)  +  M( 1/N '  +  1/N' ' ) 

J  +  K  ( N )  +  M  ( 1/N  1  +  1/N"  +  1/N'") 

J  +  K(N)  +  M(  1/N '  +  1/N"  +  1/N"’  -  3) 

J  +  K  (  N )  +  M  ( 1/N '  +  1/N"  +  2/N  "  '  ) 

J  +  K(N)  +  M( 1/N  -  4) 


=  total  number  of  carbon  atoms  in  solute 
molecule 

=  number  of  carbon  atoms  in  respective  branches 
of  branched  compounds,  counting  the  polar 
grouping;  thus,  for  t-butanol, 

N'=  N  "  =  N '  "  =  2 


TABLE  G.2 


Hydrogen  Bonding  Factors  for  Gases  in  Water  at  25  deg  C 

1 3 

(adapted  from  Hayduk  ) 


C4H10 


C2H3CL 


C3H8 


C2H6 


C2H2 


C2H4 


0.0452 


0.0031 


0.103 


2.46 


0.0326 


0.0012 


0.0359 


0.0247 


0.0053 


0.0359 


0.0071 


0.0069 


0.017  5 


0.0158 


0.0175 


0 .0134 


0.0119 


0  .0355 


0.04 10 


0.04  54 
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